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The purpose of the present paper is to study the influence of wall-echo on pressure fluctu- 
ations p', and on statistical correlations containing p' , viz redistribution 0.^ and pressure 
diffusion . We extend the usual analysis of turbulent correlations containing pressure 
fluctuations in wall-bounded DNS computations [Kim J.: J. Fluid Mech. 205 (1989) 421- 
451], separating p' not only into rapid p'^^.-) and slow p'^^^ parts [Chou P.Y.: Quart. Appl. 
Math. 3 (1945) 38-54], but further into volume (weakly inhomogeneous; P(,..2j) andp'^^.gj^) 
and surface (strongly inhomogeneous wall-echo; p'^j..^,-) and P(s u))) terms. An algorithm, 
based on a Green's function approach, is developed to compute the above splittings for 
various correlations containing pressure fluctuations (redistribution, pressure diffusion, 
velocity/pressure-gradient), in fully developed turbulent plane channel flow. The DNS 
results are then used to assess several representative second-moment closures for near- 
wall turbulence. An extended representation basis including a tensor based on the unit 
vector in the dominant direction of turbulence inhomogeneity is proposed and applied in 
representing different models in a common framework. 

Key Words: turbulence, pressure fluctuations, rapid and slow terms, wall echo, DNS, 
second-moment closure, Reynolds-stress model 



1. Introduction 

Understanding the physics of turbulent fluctuations of pre ssure p' is of major im por- 
tance, not only because of their dire ct implicatio n in noise ( Hu et a/1 12002 . 2006f) and 
excitation of immersed solid surfaces (Corcos Il964l) . but also because they appear in cor- 



relatio ns present in the transport e quations for the Reynolds-stresses and the dissipation 

is based on the Pois- 



relatioi is present m tne transport e quations tor trie Keynoids-stre 
tensor ( Pope 2000l : Jovanovid 2004 ) . Trad itionally the analysis of 



son equation for the fluctuating pressure (IChoulll945l) . which, at the incompressible flow 
limit (p ~ const = p ; Vi, a; and ^ ~ const = p, ; Vt, a;), in a nonrotating frame-of-reference, 
reads 



1 „9 , 1 „9 , , , , „ du'i, dui 
P P 
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slow 



The incomp ressible flow Poisson equation (|1.1|) suggests that ( Choul [l945l ) solenoidal 
(|Hamballl999l ) pressure fluctuations, associated with the fluctuating velocity field, are 
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generated by 2 separ ate mechanisms (|Rottalll951al f6l: iLumlevI [l97i |Piquetlll999l : iPopd 
l200d Ijovanmd3l2004l ): 

(a) the interaction of velocity fluctuations with mean- velocity- gradients called mean- 
flow/turbulence interaction terms, also termed rapid pressure fluctuations, because they 
interact immediately with an imposed mean- velocity- gradient, or linear pressure fluctu- 
ations, because the corresponding source-term is linear in velocity fluctuations, and 

(&) the turbulence/turbulence interaction, also termed slow pressure fluctuations, be- 
cause they react much slower than the rapid ones which are directly driven by mean- 
velocity-gradients, or nonlinear pressure fluctuations, because the corresponding source- 
term is quadratic in velocity fluctuations. 

This idea of distinguishing between pressure fluctuations associated with mean-flow- 
gradients and pressure fluctuations associated with turbulence/turbulence interac- 
tions only (p^x) can be applied in general to all correlations which contain the fluctuating 
press ure ( Rotta 1951al f6l: lLumlev 119781: lHaniahcl[l993 lPiauetlll999t lPopell2000l: iJovanovic 
2004. 



Choul (|1945l ) pointed out that, because of the linearity in p' of the incompressible flow 



Poisson equation (|l.ip . separate solutions can be obtained for each of the 2 source-terms , 
based on the freespace Green's function for the Poisson equation (jKatz fc Plotkinlll99l[) 
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where p{x^ t) and u{x^ t) are the pressure and velocity at point x, and the volume dt)(p, t) 
and surface d©(y, t) integrals are taken over all other points y where the pressure and 
velocity are p(f, i) and u(y, i). Notice that if (|1.2|) are multiplied by a function of this 
function can be entered into the integrals which are over y. In the case of unbounded 
flow, where 32J is very far (at infinity), only the volume integrals remain. On the other 
hand, for flow near solid boundaries, the surface integrals indicate that the unsteady 
pressure field reacts to the presence of the wall (surface integral; n is the normal-to- 
the-wall coordinate, directed outwards from the fluid volume 93). Terms related to the 
surface integra ls are usually called wall-echo terms, also known as wall blockage effect 
(|Hanialiall994[) . since for an isolated infinite plane solid boundary they can be related 
to reflection from the wall, using the r nethod of irnages (IShirlll973t iLaunder et aLlll975l : 

Gibson fc Laundeijll978t |Piquedll999l : |Pope"2000 yjovanovicll2004l). 

Starting from the seminal paper of Chou. (1945h all models ( Rottal 1951al f6l: Lumley 
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1978t IPiqued [19991 : |Popdl200d: I Jovanovid [iool for the redistribution tensor (j)^,, the 
velocity /pressure-gradient tensor Hij or the pressure transport p'u'^ appearing in the 
pressure-diffusion tensor d-^' , are traditionally composed of 4 parts corresponding to the 



(r;QJ) 



+ P(r;»)+P(s;»)+P(s;^„) 

in the the wall-vicinity is a difficult challenge ( Naka et a/J 12006). Direct 



There is at present no possibility to 
, and the simultaneous measurement 



splitting p' = p\ 
separately measure p' 
of p' and u' 

numerical simulation (dns) offers the possibility to directly compute the different terms, 
experimental uncertainty being replaced by the eventual influence of finite size of the 
c omputationa l box and of convergence of statistics. 

Kh^ (ll989^■ in the context of pseud ospectral DNS of incompressible pla ne channel 



mmi iiiijgjjii. m tne context or pseua ospectrai dns or mcomprcssiDic pla ne en 
flow teim et a/."l987': 'Mo ser et all 19991) . used a Green's functioii approach dlnce 



Courant fc Hilbert 1953: .Bender fc Orszad 1978 : Ockendon et all 2003 : Zauderei 



1926; 



2006; 



Mvint-U fc DebnathI 120071 ) to solve, as a function of y (normal-to-the-wall coordinate). 



the incompressible flow Poisson equation for p' (|1.1[) . for each parallel-to-the-wall wavenum- 
ber and The separate solution for ea ch source-term (11.11). permits the separation 
of slow pjgj and rapid p'^^,^ parts. In this way (|Kimlll989l:rChang et a/.l[l999h . the slow and 
rapid contributions to (jiij were computed ( Mansour et al. 19881 ). This procedure ( KimI 
1989t) is now used in a standard way in incompressible plane channel flow DNS, at least as 



far as p'^^y p'^^y and p'^^^ are concerned, and has also been extended to compressible flow 



studies (jFovsi et a/.ll2004( ). The term p'^^y usually called Stokes pressure ([Mansour et 



1988t IChang et a/.lll999l) . corresponds to the separately computed contribution of the 



wall-boundary condition 



corresponos 

(iPopd bood 



(11.173), p. 439) 



[fJ-9lnUn]w, n be- 



ing the normal-to-the-wall direction. Notice that this boundary-condition is associated 
with the source-term . [r/ ^l in the Poisson equation for p' , which is equal to at 



the incompressible flow limit ( Kim et al. 1987 : lMoser et al\ 19991 : Hu fc Sandhamll2001 
Hovas &: Jimenezll2006 ). This procedure was used by Chang et all ( 1999 ) to study the de- 



tailed contributions from different locations in the flowficld to wall-pressure wavenumber- 
frequency spectra. 



T-t_/J 1- —J 

evaluation of modelling proposals 


So et a/.lll991; Demuren & Sarkai1l993: Moin & MaheshI 


19981: Sauret & Vallet boOTi: Vallet 


2007). Nonetheless, to the authors' knowledge, the de- 



tailed 4-part decomposition used in the models has not been studied in detail, nor have 
DNS data been used to further split p'^^^ and p'^^.^ into volume (p^j-.trj) and p'^^.f^-^) and 
w all-echo (p'f^.^) andp'^ ^^^) terms. 



Manceau et al. 



(|200l[ ) have studied the wall-echo problem from the specific point- 

of-view of the so-called elliptic relaxation approach of DurbinI (1993), which lumps to- 
gether all wall-echo effects (both on pressure correlations and on dissipation-rate) into 
a single tensor , using 2-point correlations (in the physical space) from the dns data 



of Moser et all (199£), at wall- friction Reynolds- number Rct^ — 590. Although that 



work (Manceau et al 



20011 ) focuses on elliptic relaxation, it does provide information on 



wall-echo. The method used , although related to classical work on wall pressure spectra 
(jPanton fc Linebargeiill974l ) in that it uses 2-point correlations in physical space, intro- 
duces wall-echo by an approximate method of images, one image-channel for the upper 



f Throughout the paper x denotes the streamwise coordinate with corresponding velocity- 
component u, y denotes the normal-to-the-wall coordinate with corresponding velocity-compo- 
nent V, and 2 denotes the spanwise coordinate with corresponding velocity-component w. In !j2] 
and in the Appendix, y — at the channel centerline, while in ^and ^y — at the (lower) 
wall. 
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wall and another for the lower wall. Furthermore, iManceau et ali (|200ll ) neglect Stokes 
pressure. When analyzed under a formal Green' s funct ion framework, the approximation 
in the method of images used in anceau et ali (|200l[ ) is equivalent to using an approxi- 
mate Green's function, instead of the exact one. Based on this approach we can evaluate 
the approximation error of the method of images. 



The comparative evaluation of different turbulence closures (jPatel et aLlll985HSo et 
Il99lt IPemuren fc Sarkai 19931 ) is essential, not just for the purpose of choosing the best 
(alternatively eliminating the worse) models, but evenmore, especially in the case of 
second-moment closures, for evaluating different modelling strategies. This is particu- 
larly imgortan^Jor redistribution closures, because there exist several families of niod- 
els (I Schiestej 2007 ). and many of these, especially models including n ear-wall effects 
So et al. .1991 ), are often presented written in different tensorial bases ( Launder et 



i|ju ei aie oiieu piescmeu wiiuuen iii uuieiem uensoiiai oases uijaunuei ei an 

1975l:lGibson fc LaundeJigTSlSpeziale et a/.lll99ll:rSo et aZ.lll996l:ISugal20oi Jakirlic et al\ 



2007 ) . Contrary to the case of quasi-homogeneous turb ulence, where it is known (jSpeziale et al 



1991 



Ristorcelli et ali Il995t iGerolvmos et ali 1201061 ) that a polynomial representation 
basis of 8 tensors is sufficient for the representation of all known models, very little 
has been done for the development of representation bases in the inhomogeneous case 
(|Gerolvmos et a/.ll2004l ). 



In the present work we are interested in determining the wall-effects on correlations 
containing the fluctuating pressure. In 31] we revisit the Green's function f H2.2p approach 
( Kimlll989l) applied to DNS computations of incompressible plane channel flow ( ij2.ip . 
and we develop a procedure for separately evaluating (i j2.2.3p the volume and 
P'{s-m)^ and the wall-echo {p'l^^.^-^ and p'^^^.^-^) terms in (|1.2p . By comparison with these 
exact results we also quantify the error of the approximate method of images, as a 
function of the wavenumber f §2.2.4|) . showing that this method is a high-wavenumber 
approximation. Analytical details on the Green's functions, with particular emphasis on 
the singular problem of bulk pressure fluctuations (k = 0), are given in the Appendix. 



In S|3]we apply the algorithm ( §2.2|) for computing the p'-splitting (|1.2|) to the study of 
the contribution of different terms to correlations containing p' in low-Reynolds-number 
plane channel flow. In 21 we use these results to assess several models, representative of 
various approaches. We do not include models for pressure-st rain redistribution t^,;,- based 



on the elliptic relaxation approach, as t hese were studied inlManceau et ali ()200l[ ). but 
we study the pressure-diffusion model of Demuren et al. ( 19961 ) which uses elliptic relax- 
ation for its slow part. The tensorial representations of the differe nt 4)^ models which 
were investigated (^4.11) use different bases in the original papers (Launder et al\ 1975 



Gibson fc Laundeil 119781: ISpeziale et 0^.1119911: ISo et 0^.1119961: ISuga 



20041: iJakirlic et al 



2007^ . rendering the critical analysis of the modelling choices complicated. For compar- 
ison purposes, the various models are expressed in a common re presentation ba- 



SIS 



. -4.211, which exteiids the usual homogeneous representation bases (jRistorcelli et i 

19951 : IGerolvmos et al\\2mM\ to include an additional tensor based on the unit vector m 



the dominant direction of turbulence inhom ogeneity (§4.21) . Contrary to previous work 
( Gerolvmos et al. 2004), where a functional (Smith 197l[ ) representation basis had been 
considered, which did not allow for ex plicit projection (represe ntation) of various mod- 
els, a polynomial representation basis ( Spencer fc Rivlinlll959l ) was preferred. A priori 
assessment of different closures by comparison with DNS data is then used ft j4.3p to 
highlight strengths and weaknesses of various closure strategies. Finally, several repre- 
sentative models for pressure diffusion d^f^ (E31) are investigated by comparison with 



DNS data (g3 
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Figure 1. Comparison of present DNS-computed {Rct^ — 179; Mcl = 0.34; 193 x 129 x 169 grid; 
Tab. [1} 2-point correlations (|2.1|) containing the fluctuating pressure in wall-scaling ([i?'''^'',]"'", 



[i?<T^,]+, and [J?^Ti,]+ in the homo geneous streamwise (a;) direction; [R , ,] , [R ,^,] , and 

in the homo geneous spanwise (z) direction), as a function of distance in outer scaling 

{5~^rx and S~^rz), with reference re sults of incompressible pseudospectral ijKim et alWlQSih 
DNS computations (jMoser et a/.H 19991. J?e^„ = 178, Mol = 0). 



2. DNS computations and p'-splitting 

2.1. Plane channel flow configuration and computational method 
In the particular case of plane channel flow, the boundary surface consists of the upper 



—S and Wu at y 
and z 



+ iLy = +S, respectively), and 



and lower walls {Wi at y = —^Ly 

the periodic boundaries at x = ^-^L^ and z = zL^Lz- The methodology described in the 
present paper is independent of the particular DNS solver used. The DNS database was 
UtSv~^ = 180 (where 26 = Ly is the channel height, Ur ■= \/ fwp~^ 



generated for Rct 



is the frictio n velocity and v the kine matic viscosity) using the DNS solver described and 
validated in lGerolvmos et al. ( 2010ct ). which solves th e flow-equations in phy sical space, 
with 0(Aa;^'') discretization of the convective terms ( Gerolvmos et al. 20091 ). The flow 
is modelled by the compressible Navier-Stokes equations with air as working medium. 
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Table 1. Parameters of the dns computations [L^, Ly, {N^, Ny, N^) are the dimensions 
(number of grid-points) of the computational domain {x — homogeneous streamwise, y — nor- 
mal-to-the-wall, z = homogeneous spanwise direction); 5 is the channel halfheight; Ax'^ , Ay+, 
AycL, Az^ are the mesh-sizes in wall-units; (•)u, denotes wall and (■)cl centerline values; Ny+^-^Q 
is the number of grid points between the wall and = 10; Rbt^^ :~ u-rSu'^; u-r is the friction 
velocity; 5 is the channel halfheight; = is the kinematic viscosity at the wall; Mcl is the cen- 
terline Mach-number; Af^ is the computational time-step in wall-units; f^^g is the observation 
period in wall units over which statistics were computed; At J is the sampling time-step for the 
single-point statistics in wall-units; At is the sampling time-step for the two-point statistics 
in wall-units]. 



but for the quasi-incompressible ^ 0.34 Mach-number considered in the present 
work mean density p variations do not exceed 1.5%, and density fluctuations are neg- 
ligibly small p(.^s £ 0.25%p. For this reason the turbulent correlations obtained can 



standard incompressible pseudospectral DNS 


data (IKim et al. 


19871 Moser et al. 19991: 


del Alamo & Jimened 2003 


; del Alamo et al. 


2004; Hoyas & JimenedbOOd 20081). both 


for single point statistics ( 


Gerolymos et al. 


2010c, all SOM^jl and TOM^ appearing in 



the Reynolds-stress budgets) and for spectra of velocity fluctu ations in the homogeneous 
streamwise and spanwise directions (jGerolymos et a/.ll2010d . Figs. 12- 15, pp. 802-805) . 



This very good agreem ent with pseudospectral incompressible DNS (jKim et al. 
Hoyas fc Jimenej 2008 ) is also valid for 2-point correlations 



1987 



Ri,'h'{x,r) -.^ a'{x,t)h'{x + r,t) 



(2.1) 



containing the fluctuating pressure p', viz [7?^?^^^] + , [r'^^^'^]'^ , and (Fig.[T]) and 

fluctuating pressure spectr£0 [Epfp,]'^ and [-Epfp/]^ (Fig. [2]). 

2.2. Green's function solution of the Poisson equations 

2.2.1. ODEs for the Fourier transforms 

The solution to (|l.ip . with Neumann boundary-conditions at the walls ( Pope 2000l 
p. 439) and periodic boundary-conditions in the homogeneous directions, can only be 



obtained up to an additi ve function of time t (Ince 1926t Courant fc Hilbert 1953 
Myint-U fc Pebnathl [20071 ). We split the fluctuating pressure field as (|Kimlll989l) 



p'{x,y,z,t) = p[^-^{x,y,z,t) +p[,.-j{x,y,z,t) + p[^-^(x,y, z,t) 



(2.2) 



f second-order moments (|Gerolvmos et ani2010d . u[u'^ 

2yS[f. Fig. 10, p. 800, e^J'' := 2vd^,u[d^,u'.: Fig. 11, p. 801) 



Fig. 7, p. 797, p'u'r. Fig. 8, p. 798, 



\ third-order moments (|Gerolvmos et aL"2010j, u'^u'^u'i^: Fig. 9, p. 799) 
% Notice that the spectra given by Moscr et al. (199S) are twice the dft of Rpi^^, (in wall 
units), and have to be appropriately rescaled (|Brie:gs fc Hensonlll995l ') by (27r)~^L+ or {2n)~^Lf 
to get actual 1-D spectra E^^„. (jPopdl2000l . (6.206), p. 225), as those given bv ldel Alamo et al\ 
l|2004l ). 
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Figure 2. Comparison of present DNS-computed (J?er,„ = 179; Ma. = 0.34; 193 x 129 x 169 
grid; Tab. [l| 1-D pressure spectra in wall scaling ([£^7^/]^ and in the homogeneous 

streamwise {x) and spanwise [z) directions), as a function of wavonumbcr in waU units (k^ and 
nt), with reference results of incompressible pseudospectral f)Kim et a/... 1987. ) DNS computations 
ijMoser et a;.lll999l . Rer^ = 178, Mcl = 0). 



wher e the rapid p'f , slow p'^^^ and Stokes p'^.^-) pressure fluctuations are solutions of ( KimI 
i989HChang aijll999l) 







with boundary-conditions (|Popell200d pp. 390-392, 439-442) 

{x,y = ±^Ly,z,t) 



d_ 

dy 



P'is) 
P'(r) 


























(2.3a) 



(2.3&) 



Since ()2.3p is linear in P(„j) {m e {r, s, r}) the 3 problems can be solved independently 
to compute the 3 fields in ()2.2p . The last field in the superposition (|2.2p is introduced to 
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satisfy the wall boundary-conditions associated with the normal-to-the-wall [y) momen- 
tum equation 



dy 



y=±^Ly 



dxe 



y=±^Ly 



d v 



dy"^ 



(2.4) 



because of the no-slip wall boundary-condition u'^{x,y — ±^Ly,z,t) = 0. By p. 41) the 
field p'^^j is obviously related to the fluctuating wall-shear-strcss, and is usually called 
Stokes pressure (iMansour et aZ.lll988tlKi3ll989tlCh"ang et al. 1999) ahhoughHopi (|2000l 

p. 439) suggests the alternative term harmonic pressure, because in incompressible flow 



V'p'^ = 

The di rections x and z being homogeneous (periodic in the computational model), 
following iKim ( 1989I) . we replace the a;2;-Fourier-transformf[|of p'^^^^, Qj^) {m G {r, s, r}) 



and B' 



(t)± 

B[r)^{x,Z,t) 



+00 p-\-oc 
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+ 00 p-\-OQ 



— oo J —oo 
+00 r-\-oc 



— oo <J —oo 



(2.5a) 
{2.5b) 
(2.5c) 



in the PDE0 (|2.3p to obtain the OD^ 



5y2 



dy 



P(r) 
P'is) 
P'(r) 



(r) 







B 



(2.6a) 



(2.6&) 



(2.6c) 



The classical solution of (|2.6p by iKim using a Green's function approach, 

provides det ailed informa t ion o n the structure of the rapid and slow fields, and was 
extended bv lChang et al. ( 1999t ) to study the detailed influence of the sources, both in 
type (slow and rapid) and y-location in the channel. In the present work we extend 
this analysis by identifying and studying the wall-echo influence. Details on the Green's 

functions, and discussion of the singular case k — Kx — — are found in the 

Appendix. 



f This is however no longer true in the compressible flow case l|Fovsi et aLll2004l l. 
I The Fourier-transforms p[ m): Q[m)j B'(^m)+ £ C are computed using standard dft (discrete 
Fourier transform) techniques (jBriggs fc Hensonlll995l) in the periodic directions x and z, with 
maximum computable wavenumbers k^;,^^^ = 7r(Aa;)~ and Kz^ax = tt{Az)~^. 

% PDE: partial difi'erential equation; ODE: ordinary diflterential equation 
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2.2.2. Exact solution of 

The exact soluti on to (I2.6p is fiiven in [Kiml(ll989l) . and has been widely used (^Mansour et al 
1988HChang a^.i.lQQQ; .Fovsi et a/.ll2004l) . It reads f COTTl COl]) 



P'ir) 
P'is) 



{f^xj f^z, t) 



1 r 

- 2 -tJy 



2 « 



GKiniiy.Y; k) 



'(r) 



(2.7) 



where the Green's function GKim{y, Y; k) (Tab. [5]) is defined by (jA 961IX34)) . and satisfies 
the correct homogeneous Neumann boundary-conditions on both walls (Fig. [3]), and 
K := \/ kI + e R>o (|2.6cp . The Stokes pressure is given by 



P[r){'^x,y, Hz, t) = (7Bc(y; «, B'(^r)- (l^x, Kz,t), (k^, Kz,t)) 



(2.8) 



where the function qBciy, k, B-,B+) (Tab. [2]) which solves for the harmonic field p'^^j is 
defined by (jA 2ip . As discussed in the Appendix ( ijA.2p . the case k = is singular, and 
a solution exists iff the compatibility conditions 



Q', 



-xyz 



_1 r 

2 y 



{k^ ^ 0,Y,K, ^ 0,0 dY 



(2.9a) 



B[^-^_{K^ = 0,K,^0,t) = B[^^_{K,^0,K, = 0,t) ^ B[^^_"\t) = B[^^_ 



'it) 
(2.96) 



where (•) denotes the volume averag el (bulk average) and (•) is the surface averag^ 
in the homogeneous directions x and z. If the co mputational box is large enough in the 
homogeneous directions for the ergodic hypothesis (jMonin fc Yag lcm]ll97ll p. 243-256) to 
hold, ()2.9p are satisfied, at least approximately. In practice, the compatibility condition 
(|2.9ap was enforced when computing the Fourier-transforms (|2.56l) . The form (jASip 
of qBc{y]K, — 0, i?_,i?+) (Tab. [2]) automatically enforces the compatibility condition, 
by taking the average gradient \{B^ + -B+). In practice B'^^^_{kx = 0, — 0,t) ~ 

B'^^s^^{kx = 0, K2 = 0,t) ~ in (j2.9&p . especially with increasing box size, so that taking 
the average is a good choice. The additive con stant Qbcq up to which (jl.ip with Neumann 
boundary conditions at the walls (|Ponell20nnl p. 439) can be solved was chosen to satisfy 
the constraint p'^{t) := p'^^^ (t) = 0, (ie constant bulk pressure) and was included in the 
Stokes field p'^^^ . 



2.2.3. Volume and wall-echo terms 



The main purpose of the present algorithm is to distinguish in (II. 2p the contributions 
of the volume integral (volume terms denoted by {■)<xi) from the contributions of the 
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G^^Av, y; . / 0) = „ coshKL,^|y^y|)l + coshK, + y)] ^ 



G±[y,Y-K^O) = G<r,{y,Y-K^Q) 

I _Lc.u u\ 5+cosh[K(iLy + y)]-B_cosh[K(iLj,-y)] 
qBc{y;tij^O,B-,B+)^ . (jAllJl 



f+i^y , , 

compatibility condition: / Q[Y; k = 0)dY = (|A 24|l 

GKim(y,y;«; = 0) = G±(y,Y;«; = 0) = G<£,(y,F;A. = 0) = jXH 

(7Bc(y;«: = 0,B_,B+) = iB.y + + Qbco 

gBc±(y;«; = 0,B±) — (IX36|l 

Approximate method of images f i]A.3P 

GMWL(y, y; «: / 0) =Gv{y: F; k / 0) + Gu._ (y, Y; k / 0) + G^^ {y, Y;k^O) (IXITBi 

/ r \ e-''^y{l + e-''^y) cosh Kly j,-^-^ 

£mwl(kL;,) ;= )— — ^-r^ (|A 44|l 



Table 2. Green's functions for the solution of the generic modified Helmholtz equation (|2.6ap . 
q"{y;K) — <l{y',n) = Q{y',K) Vy € {—^Ly,+^Ly), with boundary-conditions (Fig. [3| 
q'{y = ±iLy;K) = B±{k,) (|2.66|) . and for the two halfspace problems (Fig. [3]) defining the 
wall-corrections of the method of images (Fig. [Sj; notice that the solution for k = is obtained 
up to an additive constant ^Bro (|A 31|) . 



surface integral (wall-echo terms denoted by {■)w), ie to extend (|2.2I) as 

" V ' 

+ p[s-^){x-,V, z,t) + p\^.^){x,y, z,t) +p\^^{x,y,z,t) (2.10) 

v ' 
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Figure 3. Instantaneous pressure fluctuation sources (sum [Q(r) + Q'(s)]^ '^f rapid and slow 
sources in waU units) at x = (dnS; Re-r^ = 179; Mcl = 0.34; 193 x 129 x 169 g rid; Tab. [l]), 
and boundary cond itions for the exact problem (|2.6p whose solution (|Kiml [1989 ) is given by 
GKim(j/, Y; k) ()A 96|) . and for the two virtual halfspace problems with only one wall present defln- 
ing G±{y, Y; k) (|A 17|1 . which deflne the local wall-corrections appearing in the high-wavenumber 
approximation of the method of images (|Manceau et aZ.l200ll '). corresponding to the approximate 
Green's function GMvi-L{y,Y; k) (|A 41c[) . 



where p'^^.t^^ corresponds to the contribution of the volume integral in (ll.2ap . p'^j..^-, cor- 
responds to the contribution of the surface integral in (|1.2ap . p^g.sjj-j corresponds to the 
contribution of the volume integral in (|1.2&[) , and p^g.^-j corresponds to the contribution of 
the surface integral in (|1.2&p . From a modelling point-of-view, the volume terms Etnd 
^'(s-aj) correspond to a quasi-homogeneous (weakly inhomogeneous) part, while the sur 



face terms p^^.^,) and p'^^^-w) ^^'^ strongly inhomogeneous wall-blockage effect (jHanjalic 
11994 . 

Working directly with the surface integrals in (|1.2p would have been a complex task, be- 
cause they are implicit, ie they contain the value of the corresponding field at the wall. On 
the other hand, it is straightforward to evaluate the volume integrals in (II. 2[) . which con- 
tain only the sources and the freespace Green's function G<x3{y, Y; k) which solves (|2.6ap 

replacing the boundary-conditions (j2.6&P by the requirement that hm|j,| >ooP(r tjj) ^ 

lim|j,| s.ooP(i..2j) = 0; the solution of the hypothetical problem where p' is generated 

by a distribution of sources in y e [—^Ly, +^Ly] with the walls absent. Straightforward 
computation fi jA.1.31 i]A.2.3l) gives the freespace Green's function G<s{y,Y; k) (jA 13[ 
lA 34|) for the 1-D case (Tab. ^ which satisfies 



and since by definition ([2T0|) 



+ 2 



P(r;w) 
Pis;w) 



{I'^x, y, i^z, t) 



P(r) 
p'is) 



G<n{y,Y;K) 



{K^,y,K^,t) 



i^^xi Y, Kj., t) dY 



(2.11) 



^(r;^) 
P{s;<XS) 



{n^,y,n,,t) (2.12) 



we may readily define, using the Green's functions (Tab. [5]) GK-,,y,iv.Y: K) (jA9M lA 34p 
and G<s{y,Y\K) (|A 13[ 1X341) . the wall-echo Green's function 



Gw{y,Y]K) :=GKimiy,Y;K) - G^{y,Y]K) 



(2.13) 
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corresponding to the surface integrals (wall-echo) in (|1.2 



P(s-w) 



_1 r 
2 » 



(2.14) 



Notice (Tab.H]) that since G]^i^{y, F; k = 0) = G'aj(?;, F; k = 0) (|A34|) we have G„(?;, F; k ^ 
0) = 0, ie wall-echo applies only on fields varying in at least one of the homogeneous 
directions x and z. 



2.2.4. The approximate method of images 

With the above developments we can directly evaluate the wall-echo terms for the 
slow and rapid wall-echo terms in each of the correlations coi itaining p' (jlSl). We st udy 
the method of images i iot only because it has been used bv iManceau et al. ( 200l[ ) to 
analyse the D NS data of Moser et al. i 19991 ). in the context of elliptic relaxation models 
( Durbinl[l993h . )ut al so because it is often used as a closure-building tool in near- wall 
turbulence modelling (|Shiilll973tlLaunder et aLlll975tlGibson fc Laundeiill978l ). By com- 
paring with the exact solution obtained in the present work, it is possible to evaluate the 
approximation error of the method of i mages. 



The approximate method of images (jManceau et ali [20011) simply adds two mirror 
images of the channel, one above the upper wall and another below the lower wall (Fig.|3]), 
to account for the presence of the walls, postulating that an approximation to the solution 
is obtained by using the freespace Green's function G<2j{y,Y; k) in the extended domain 



y G [- 



3 r 

' 2 y " 



\Ly]. Taking into account the mirror symmetry of the ghost channels with 



respect to the corresponding wall (Fig.[3|), this is equivalent f ^A.3|) to using an appropriate 
Green's function Guwhiv, Y] k) (Tab.[2|) on the actual sources between the channel walls, 
ie for y £ [- 



±r -i-i 7^ 1 

2 V 2 V\ 



-P(r;MWL)J 
P(s;mwl)] 



By (|A41e|) 



+ 2^y 

1 r 

2 y 



GMWLiy,y; n) 



{K,,Y,n,,t)dY (2.15) 



Guwdy, Y; k) = G<i}{y, Y; k) + {y, Y; k) + G^^ (y, Y; k) 



(2.16) 



where Gw± {y, Y; k) are the wall-echo Green's functions corresponding to the virtual 
halfspace problems with only one of the walls present (Fig.[3|). Notice that this is also th e 
method used in standard wall-echo Reynolds-stress modelling (jGibson fc Laundeilll978l ). 
and (Tab. [5]) that the halfspace problems are solved exactly f ^A.1.4[ ^A.2.4[) . Inherently, 
this is tantamount to assuming that there is no interaction between the echo effects of 
the upper and lower walls, since in the case of an isolated wa ll (halfspace problem; Fig.[3l) 
the method of images yields the exact solution (|Ponell200nl pp. 439-442). 

The approximation error of the method of images comes from the fact that the direct 
influence of the upper wall, approximated by the halfspace problem Gw+ {y, Y; k) (Fig. [3]), 
induces a nonzero gradient dy[p'^j.._^yp'^^._^-^Y{Kx,y = —■^Ly,Kz,t) ^ at the lower wall, 
and that there is no feedback from the lower wall to correct this (and vice versa). The 
nondimensional appro ximat i on err or can be evaluated ( ^A.3[) by comparison with the 
exact solution (jA 9bp of Kim ( 1989( ) , and is a function of the nondimensional wavenumber 



KLy = 2kS (|A44|) . For large nondimensional wavenumbers uLy, ie structures with small 
streamwise and spanwise extent compared to the channel height, this error is small, and 
rapidly decreases with increasing wavenumber (Fig. |4|). 

Let £k. '■= 27rK~^ be the representative size (wavelength) in the xz plane of structures 
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Figure 4. Comparison of the exact Green's function (|A9fe|l . G\^^m(y ,Y : k) l|Kiml Il989f) . 
with the method-of-images approximation (IA41|I . Gt.iKi,{y,Y;K) ijManceau et al\ |20 0ll). for 
Y £ { — T^Ly, — -jLy, 0} and uLy G {1,|,2, 3}, plotted as a function of yLy^, and log-plot 
of the upper bound estimate of the error made by the method-of-images approximation (|A 44|l 
as a function of the nondimensional wavenumber uLy — 2kS {Ly = 25 is the channel height) . 



corresponding to k := y/ + p.6cp . For nLy ^ 3 ~ iS the approximation 

of GKim{y,Y; k) by GMWh{y,Y; k) (Tab. [2|) is satisfactory (Fig. |4]). For the present case 
with Rcr^ ~ 180 this corresponds roughly to ^ which includes the wavenumbers 
where most of the energy of the p'-spectra is contained (Fig. [2]) , implying that the method 
of images is a satisfactory engineering approximation. However, for larger structures, eg 
KLy ~ I <^=^ £k. ~ 85, such as the superstructures observed in higher- i?eT-„ wall 
turbulence ( Hutchins fc Marusic 200?! : iBalakumar fc AdrianI 2007 ). the approximation 
error is roughly 10%, and then grows exponentially with increasing size. 



2.2.5. Interaction of wall-echo between walls 

Since the freespace Green's function Gaj(y,i^; n) ()A 13| lA34|) is the same for all prob- 
lems (Fig.[3|), the differences (Fig.|4]) between the exact solution G'Kim(2/, Y] k) (|A96[) and 
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Figure 5. Comparison of the exact wall-echo Green's function Gw{y = —\Ly, Y; k) (|2.16p . for 
the computation of the solution at the lower wall y = —^Ly, with the wall-echo Green's function 
Gw_ {y = —^Ly, Y; k) (|A 411) . of the corresponding halfspace problem with only the lower wall 
present (Fig. [3]), for KLy £ {1, 2, 3}, plotted as a function of YLy^ , and log-plot of the ratio at 
the wall (y — Y — —^Ly) oi the exact wall-echo Green's function Gw{y ~ —^Ly,Y = —^Ly;K.) 
1)2. 16[) on the corresponding freespace Green's function Gnjiy ~ ~\Ly,Y — —^Ly;K.) ([A 13)) 
plotted against the nondimensional wavenumber nLy = 2k(5 {Ly — 2S is the channel height). 



the method-of-images approximation GMWh{y,Y;K) ()A4ip . correspond to increasingly 
strong interaction between the two walls with decreasing nondimensional wavenumber 
KLy = 2kS, leading to amplification of the echo effect. To further explain this phenomenon 
notice the identity 



i lXTal i 



G^iy = -lLy,Y;K^O) ^ 

lA 17a^ 



2k 



G-{y 



\Ly,Y; K ^ 0) - Gaj(y = -\Ly,Y- n^Q) 
^^G^_{y^~\Ly,Y-K^Q) (2.17) 



ilAlT 



which implies by ()A 3)) that for the halfspace problems ( ^A. 1.4)1 the wall-echo, at the wall 
[y — —^Ly), exactly equals the volume term at the wall. 



{x,y 



1 J- ,^ |2l3jA3 , 



2 : ^7 



(2.18a) 



_)ix,y ^ ~\Ly,z,t) ^^^^^p'f^^,^^{x,y = -\Ly,z,t) (2.186) 



This result 02.18)1 following from the equality of the corresponding Green's functions 
()2.17p is valid independently of the particular sources Q{y\K) ()2.6ap . 

At high nondimensional wavenumbers K,Ly the exact wall-echo Green's function G^ {y 
— ^Ly,Y;K) ()2.16)) . for the computation of at the lower wall {y = 
mately equal {uLy — 3; Fig. [5]) to the wall-echo Green's function Gw_, {y 
()A 4ip of the halfspace problem with only the lower wall present (Fig. [3]), except at the 
upper part of the channel Y G [jLy, \Ly\, where the influence of the upper wall on the 
Green's function used for the computation of at the lower wall {y = — \Ly) is felt. As 



-iLj,), IS approxi- 
2 ' 
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ife^„, A/cL iV,. xNyX L, L„ A.t+ Ay+ Ny+^i„ Ai/+ Az+ At+ At+_, 

179 0.34 193 X 129 x 169 ijrS 2S ijrS 11.7 0.22 20 4.9 4.4 6.47 x 10"^ 1648 6.47 x IQ-^ 




Figure 6. Instantaneous fluctuating pressure in wall units [p']"*" obtained dire ctly (Rct,,, = 179; 
McL = 0.34; 193 x 129 x 169 grid; Tab. [T| by the compressible DNS solver (jGerolvmos et al.\ 
l2010d ) compared to the superposition of the three fields, rapid [pj^)]''' and slow [p'{s)]^ obtained 
by the Green's function solution (|2.7p . and Stokes field [p'(t-)]''^ (|2.8p . in wall units (50 contours 
in the range [—5, +5] on the lower wall, the outflow a;-periodic interface, and the z-periodicity 
interface) . 



the nondimensional wavenumber uLy further decreases {ie the corresponding structure 
size £k5~^ :— 2ttk'^S~^ increases) the influence of the upper wall is felt further down 
(Fig. [3]). Of course there is an analogous effect concerning the influence of the lower 
wall on the upper wall. As the nondimensional wavenumber KLy decreases this mutual 
interaction between the two walls amplifies the echo effect, the amplification growing 
exponentially with structure size (Fig. [3]). 



3. Analysis of pressure correlations 

The algorithm ( ii2.2.3p for separating the rapid and slow contributions to p' ( H2.2.'2\\ into 
weakly inhomogeneous (p'^^.^j) and ^(s.gj')) and strongly inhomogeneous wall-echo terms 
(p[t-w) ^^'^ P[s-w)^' applied to a low- Reynolds- number [Rcr^ ~ 180) well-resolved 
DNS (grid 193 x 129 x 169; Tab.[T]). The simulation was started at t = io by interpolation 
of a well- converged simulation on a coarser grid (grid 129 x 129 x 129; Tab. [T]), and 
continued for ~ ^ct — 681. During this interval i+ G [io^j^GFtJ statistics for Ui and 
u'ju'^, necessary for the computation of the source-terms in (I2.3ap . were computed with 
sampling at every iteration (At+ = A<+ ~ 6.47 x 10"^ /+ ~ 154). Then the 

Green's function algorithm for p'-splitting (|2.10P was applied, and statistics of pressure 
correlations were computed for an observation time tjgg = 1648, with sampling every 
100 iterations (At+^ = 100Ai+ « 6.47 x IQ-^ ^ f+J'k 1.54). 

3.1. Fluctuating Pressure Field 
The instantaneous fluctuating pressure p'-field reconstructed by the Green's function 



approach (|2.71 12. 8p according to the splitting (I2.2P agrees quite well (Fig. [S]) with the 
instantaneous p'-field directl y computed (without m aking use of Green's functions) by 
the compressible DNS solver (jGerolvmos et aLll2010cl ). 

Both the directly computed p[^^ and the rms values ([P(^)]rms, [P(s)]rms, [p'(T-)]rms) of 
the p'-splitting (|2.2p compar e satisfactorily (Fig. [71) with standard re sults from incom- 
pressible DNS computations ( Kim et al. 1987; Hovas fc Jimened 2008 ), both in the wall 
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Figure 7. Comparison of rms-values Prms := V p'"^ of the fluctuating pressure p and of 
the fields p'^^j + p'^^j + p'^^j = p (|2.2|) from the present DNS computations (Rct^ — 179; 
McL = 0.34; 193 x 12 9 x 169 grid; Tab. [1]) with reference results of incompressible pseudospec- 
tr al (iKim et al\ 1987 ) dns computations of Hoyas & Jimenez (2008, Rer^ = 186, Mcl = 0) and 
of|M oser et al] (|1999| . Rgt^ = 178, Mcl =0), in wall units, plotted against the nondimensional 
distance-from-the-wall in inner (y'^) and outer {5~^y) scaling. 



and the outer regions (Fig. [7]). There is, nonetheless, a difference in the level of Stokes 
pressure [p'^^^Jmis near the wall (Fig. [T]), which should be further investigated. However, 
this discrepancy has no impact on the a priori evaluation of various closures O, since 
all models ( H4.ll S31) neglect the small contributions of the Stokes pressure p'(^^y al- 
though, as discussed in i j4.3l the contribution of Stokes pressure to 4>xy is dominant in 
the near- wall region (y+ ^ 10). 

The new results (Fig. [5]) in the present work concern the further splitting [p' — P^^.ijj) + 
P[-[-w) ~^ P[s-v) ~^ P'(s-w) ~^ p[t)) volume and wall-echo terms (I2.10p . The rms-levels of 
wall-echo ([p^^.^^lrms, [P(s;„)]rms), at the wall (y = 0; Fig. E]), are approximately equal to 
the corresponding volume terms ([P(r.(i;)]rms, [^(s-aj)]''™^)- This implies f ^2.2.5p that the 
dominant contributions to the spectra (j2.56p of the source-terms Q'f^^'^R{Kx,y, Kz,t) and 
Q(s)('*^' '^^^ ^) occur at sufficiently high nondimensional wavenumbers nLy = 2k,5 

(Fig. [5]), for the interaction between the wall-echo from the upper wall and the wall-echo 
from the lower wall to be negligibly small. This approximate equality of wall-echo and 
volume terms holds up to y"*" ^ 2 (Fig. [8]). Further away from the wall, both wall- 
echo terms ([P(i..^)]rms, [P(s ,i,)]rms) dccay with increasing distance from the wall, always 
remaining much higher than the Stokes pressure [p'^^jjrms (Fig. [8]). In the buffer and outer 

regions (y ^ T§o'^ '^^^ ^ known predominance of the slow volume 

term [p'^^ 

.Qjj]rms over the rapid volume term [p|j..Qjj]rms 

(|Kimlll989l: IChang et aZ.lll999t) is 
evident (Fig. [5]). Nonetheless, it is essential, from the point-of-view of near- wall modelling. 
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Figure 8. Rms- values of the 5 fluctuating pressure fields in the decomposition 
P' = P'{r;?)) + P'(r;ii.) + p[a;m) + p[a;w) + p'(t) l|2-10p from the present DNS computations 
(i?er„ = 179; Mcl = 0.34; 193 x 129 x 169 grid; Tab.[T]), plotted against the nondimensional dis- 
tanco-from-thc-wall in i nner (j/"*") a nd outer (S^^y) scaling, and p'^^s from various DNS databases 
(Moser et a/.., 1999. : .Hovas fc Jimenez., 2008. : .Gerolvmos et a/...2010G ') . 



to notice that for y+ £ 10, b(,.2j)]rms ~ [P[s-<XS)\™^ - \P'(r;w)\™^ - b'(s;«,)]rms > [P(.^)]rms 

(Fig. IHl, ie for consistent near- wall modelling all of the 4 terms must be accurately 
modelled. 

The above quantitative results (Fig. [5]) are also observed qualitatively in the instanta- 
neous levels (Fig. 19]) of the 5 terms in the p'-splitting (|2.10l) . Comparison of instantaneous 
levels of the slow [p'^g-^] and rapid (P(r)) terms (Fig. [5]) clearly shows that both mechanisms 
of generation of p' (I2.3al) are of the same magnit ude at the wall (lower wall; Fig.lH]), while 
p'^^-^ is the main mechanism ( Chang et a/.lll999l ) further away from the wall {z = const 
plane; Fig. E]). The same observations apply to the corresponding volume terms, p'^^.s^j-^j 
and (Fig-IHl)- The wall-echo terms, p'^^.^^ and p'f^^.^y at the wall (lower wall; Fig-IH]), 

are approximately equal to the corresponding volume terms, p'f^j...^-^ and P^g.^j-j , but rapidly 
decay away from the wall (z = const plane; Fig. [S]). Finally, the Stokes pressure p'^^-^ is 
substantially lower than the other terms and rapidly decays away from the wall (Fig. IH]). 

The results (Figs. EHll) on the p'-fields (|2.10p provide guidance on the relative im- 
portance of each of the 5 p'-fields. Nonetheless, the transposition of these results to the 
decomposition of correlations containing p' is not always straightforward, especially in the 
near- wall region (y+ ^ 10). For this reason we study in detail f ii3.31 WiA[ i j3.5p the cor- 
relations containing p' which appear in the transport equations for the Reynolds-stresses 
(». 

3.2. Reynolds- stress transport 
The equations governi ng the R eynolds -stress tenso r are central in single-point closure 



turbulence modelling ( Lumlevlil97& .Hanialic 1994). The exact transport equations for 
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iJe^„, A/cL iV.,. X N.y X L,. L„ A.t+ Ay+ Ny+^w Ay+_ Az+ At+ At+_ 

179 0.34 193 X 129 x 169 4m5 2S ijrS 11.7 0.22 20 4.9 4.4 6.47 x 10"^ 1648 6.47 x IQ-^ 




Figure 9. Instantaneous fluctuating pressure flelds (|2.10l) from DNS computations {Re-r^, = 179; 
McL = 0.34; 193 x 129 x 169 grid; Tab.©, Stokes b'(^)]+ ^M, and rapid and slow [p'(^)] + 

l|2.7|l . and their decomposition into weakly inhomogeneous volume fields, [P(i..(jj)]^ and [pjg.tjj)]^ 
l|2.1H) . and strongly inhomogeneous wall-echo flelds, [p'j^. „)]"'' and [pj^.^,)]^ (|2.14p . in wall units 
(50 contours in the range [—5, +5] on the lower wall, the outflow x-periodic interface, and the 
2-periodicity interface) . 
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the Reynolds-stresses in incompressible flow read (jPope 1200(1 pp. 315-320) 



convection diffusion + d*^' + 

'■J '■J ''J 



du' du'A ( 



redistribution 4>ij production 



dissipation pe.^J-' := 2p.dx(U^dxiu'^ + piidxjU'fdx^i + dx^u'^d^^u'-) 

Introducing an exact term for the viscous diffusion of the Reynolds-stresses {d\^'^ ) instead 

of the original term ('^ij'*), this exact equation is recast in the following mathematically 
equivalent form, which is actually modelled (|PoDell2000l (7.191), p. 319) 



dpu'-u'- d _ d du[u'- 



convection C,, diffusion d^") -I- d'f? -f d^^^ 



redistribution 0y production Pij 



d du':u', dr' dr'n \ 



dissipation pe^^'' := 2p,dxeu[dxiu'^ 

Convection Cij, production Pij and viscous diffusion are exact terms, while all 

the other terms inj, and e^^^) require modelling. Pressure-diffusion and 

redistribution (f)ij can be grouped together into the velocity /pressure-gradient correlation 
tensor 11^ 



The velocity/pressure-gradie nt tensor 11 ,;,, which is the term appearing in most of the 
original early developments (|Choi]lll945h . tends to as y — >• (in the viscous sublayer, 
where ci|^^ and ipij cancel one another). In homogeneous turbulence (ht) dff^ = <^=4> 

0^.i — 11,1 . 
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Figure 10. Distributions of the 5 terms in the incompressible p'-splitting (|2.10|) of pressure 
transport, p'^^.^jii^, p^^.^^u'^, p[^.<3j^u'-, p[^.^^u[, and p'^^^^u'-, from the present DNS computations 
(i?er„ = 179; Mcl = 0.34; 193 x 129 x 169 grid; Tab. [T]), distributions of pressure transport 
p'u'f from various DNS databases (jMoser et a/.lll999l : iHovas &: Jimened [20081 : iGerolvmos et al\ 
l2010d l. and corresponding contributions to the pressure-diffusion tensor d[^\ in wall units, 
plotted against the nondimensional distance from the wall j/"*" (in fully developed incompressible 
plane channel flow d'i'J = d^J = d'fj = d'fj = 0). 



3.3. Pressure transport p'u[ and pressure diffusion d^^ 

The p'-splitting (|2.10p applied to pressure transport p'w ■ and pressure diffusion d-^^ p. 11 
13. 2p indicates (Fig. [TUl) that the slow volume term p'^^.^-^v' is the principal contribution 

to the normal-to-the-wall transport p'v' , and evenmore to d^yy (for fully developed plane 
channel flow d^y '^^^F^ —2dyp'v'), except in the near-wall region {y~^ S 10) where 
all terms are of comparable importance (Fig. [T0|) . Concerning d^xy (for fully developed 
plane channel flow d^xy ''^^^P^ again the slow volume term d^J'y^''^-' is the most 
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Figure 11. Distributions of the 5 terms in the incompressible p'-sphtting (|2.10|) of pressure-s- 
train redistribution, 4'['j''^\ <t>i'j'™\ a-nd ^'j', from the present DNS computations 
{Rbt^ = 179; McL = 0.34; 193 x 129 x 169 grid; Tab, [p, an d distribu tions of pressure-s- 
train redistribution from various DNS databases f(Moser et al\ \T99Q: Hov as fc Jimene3l2008l : 
iGerolvmos et a/1l2010d ). in wall units, plotted against the nondimensional distance from the 
wall (in statistically 2-D plane flow (jj^x = 0t/z = 0). 



important contribution in the buffer and outer regions, whereas in the near-wall region 
{y~^ ~ 10) all terms are of comparable importance (Fig. [10]). 

For modelling purposes it is important to notice that, although the slow volume term 
^(p,s,2j) ^j^^ main contribution (and a reasonable approximation) to d[j^ (Fig. [TO]), in 
the buffer and outer regions (y"*" ^ 10), this is not the case for pressure transport p'u'^, 
especially for p'u' (Fig. [T0|). where, for 5 ^ j/+ ^ 30, ^(s.gj)'"' is of the opposite sign with 
respect to p'u'. 

These remarks imply that a consistent model for pressure diffusion d[^^ , in the buffer 
and outer regions (y^ ^ 10), can be built using only slow volume terms. However, such a 
model for d-^' (y+ ^ 10), based on d\j''^''^^ only, cannot correspond to (or be built from) 
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a satisfactory model for pressure transport p'u[ (Fig. [TU)) . The near- wall modelling of 
pressure diffusion (f^^ (y"*" ^ 10) is more complex, since (Fig. [Til it must contain rapid 
terms and wall-echo terms (Fig.[TU]). To the authors' knowledge|Vallct 2003) no closure 
for d!f^ satisfies all of these requirements (iDonaldson 1969 : Hir^ 19691: iDalv fc Harlow 



1970l:TLumlevlll978l: lF^ ll993t iDemuren et allliggeTlSauret fc Valletll2007t IValletll2007l) 



3.4. Pressure- strain redistribution 
The application of the p'-splitting p.lO|) to pressure-strain redistribution 4)ij p.l) 13. 2p 

^f°^ and 4f 



shows again that the volume terms (0^^'^'' and 0^^'^'') are the dominant contribution to 



(pij in the buffer and outer regions (y+ ^ 15; Fig. [TT|) . Notice that ip'ti^^ > ^iz'^'' for 
2/^ ^ 15 (Fig. [TT|) . contrary to the other components for which the slow volume terms 
are the more important contributions to (jjij in the buffer and outer regions {y'^ ^ 15; 

Fig.in]). 

In the near- wall region (y+ ^ 10) the wall-echo terms {(t>[j'^'' and 4>l^j''"^) are of the 

same order-of-magnitude as the corresponding volume terms {(p^'j''^^ and 4>'fj''^\ respec- 
tively), satisfying an approximate equality at the wall, which (of ^2.2.5p implies that the 
energy-containing nondimensional wavenumbers of the source-terms (I2.3a[ 12. 5M I2.6ap 
are sufficiently high for the interaction between upper and lower wall to be negligible 
(Fig. [5]). Furthermore, for the low- Reynolds-number case studied in the present work 
(i?e-r„ ~ 180) the Stokes pressure term ^^J-* is of the same order-of-magnitude, in the 
near- wall region (y+ ^ 10; Fig. 1111) . as the other terms of the p'-splitting (|2.10p . and 
is particularly important for the shear component (j)xy (Fig. [TT|) . At the wall (y+ = 0; 

Fig. [Tl]) [(/)i^'']y=o accounts for ~ 25% of [(j)xy]y=a- For a plane wall _L Cy, at the wall 
(y = 0), [4>a:y]y=o ''^^F^ \p'dyu']y=o (because v'y^Q ^ =^ [<9a;W%=o = 0). Therefore, 
the high level of \<i)"xy\y=Q relative to the other terms of the p'-splitting (|2.10l) implies that 
at the wall the fluctuating Stokes pressure [p'^^-|]y=o is well correlated with the fluctuating 
wall-shear-stress [t^j,]j,=o = Ai[i9i/w']i/=o- Notice that none of known models for 0y ( i)4.ip 
gives the correct viscous sublayer behaviour of <l)xy {y'^ ^ 5; Fig. [TT|) . 

3.5. Velocity/pressure-gradient Ily 
The velocity/pressure-gradient correlation Ilij (|3.2p is exactly the sum of pressure diffu- 
sion d^j' fi j3.3p and pressure-strain redistribution 0^ ft ^3.4p . By definition (|3.2I) [Hy-Jy^o 
because of the no-slip wall boundary-condition for the fluctuating velocity [u^]j,=o — 0, 
implying [4>ij]y^o = — [dj-^'']j,=Oi so that very near the wall we expect that Hij := (j^ij+d'f^^ 
will be asymptotically approaching 0. Observation of DNS data (y"*" ^ 5; Fig. [T^ indi- 
cates that the rate at which XI^^ approaches as j/+ is much slower compared to 
the other components {Jlxxi ^xy^ and H^j,). 

In fully developed incompressible plane channel flow d^ix = d^zz = by (|3.2p , implying 

that Tixx — (pxx and Hzz — (j^zz, but also that Hu ~ d^yy (since by the 

incompressible fluctuating continuity equation dx^Ug — 0). For the three components, 
^xx, ^xy and Tlyy, the dominant contribution from the mechanisms of the p'-splitting 
(|2.10P comes from the slow volume terms {Tl^xx^\ niy^\ and Ilyy^'*), the rapid volume 
terms (Ili'i^'^^ , Iliy^^ , and Hyj)^^ ) being the main remaining contribution (Fig. [12]) . On 
the contrary, all of the 5 mechanisms of p'-generation in (|2.10p contribute, with the same 
order-of-magnitude in the near-wall region (y+ ^ 15; Fig. [T2|) to Hzz (= (t>zz in fully 
developed incompressible plane channel flow), this difference between IIz^ and the other 
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Figure 12. Distributions of the 5 terms in the incompressible p'-sphtting (|2.10|) of veloc- 
ity/pressure-gradient correlation, n'j'^\ n'j'™', 11 



(s;93) 



n 



(s;uj) 



and , from the present 



DNS computations {Rer^ = 179; Mol = 0.34; 193 x 129 x 169 grid; Tab. [T]), and distribu- 
tions of velocity/pressure-gradient correlation 11^^ f rom various DNS databases (jMoser et al\ 
119991 : iHovas fc Jimenezll2008l : iGerolv mos et ffllll2010d ). in wall units, plotted against the nondi- 
mensional distance from the wall j/"*" (in fully developed incompressible plane channel flow 



llza; — — 0, rixa; — 



and lizz = ' 



components being related to the slower rate at which Il^z goes to as j/"*" — !• (y"'" ^ 5; 

Fig. [ED. 



4. A priori evaluation of SMCs for the pressure correlations 

Pressure correlations p.2p appearing in the Reyno lds-stress transport equations (|3.ip 
are of fundamental importance in SMCs (ILaunder et al . 197^; Hanjalic 1994). Evaluation 
in complex flows around co mplex geometri es (jGerolvmos et al .I I2OO2L 120103) ■ confirms 
the generally accepted view (iHanialic that the rapid redistribution terms control 

separation fiVallet.2008h . but also suggest that pressure diffusion may be influential not 
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only on the entrainment at the edge of shear layers ( Demure n et aL* [l994l : Suga 2004h 
but also in controling reattachment in the near-wall region " (|vallet,2007 l). 

In the present work we assess a priori SMCs of the pressure correlations used in 
RSMs, separately evaluating the closure for the quasi-homogeneous volume terms and 
the strongly inhomogeneous wall-echo terms, following from the splitting (|2.10p of p' 



p'u[ 



1 

2^ 



1 

4-77 



,jj dik d^e If - x\ 



L fff d^u',u',dx>{f,t) 



[p 



[p/u'l(r;™) 



2j 'dudTk |y-a;| 



[p'u 



(4.1a) 



(4.16) 



where the volume integrals in (|4.ip are obtained by multiplying (|1.2p by 2S^j (|4.1ap or 
by u- (|4.16p . and ensemble-averaging. 

Every order-2 tensor A e E'^^'^ has 3 invariants (lRivlinlll955l l. 



Ia 
IIa 
IIIa 



trA 

i(tr2A - trA^) = 1 {{Au)^ ~ AtmAe^) 
dct A = eijkAiiA2jA^k — eijkAiiAj2Ak3 

) 



=1 (Ia3 + SIaIIa 



2trA-^ - 3trA trA^ + tr^A 



trA^ =IX - 2IIa 



(4.2a) 
(4.26) 



(4.2c) 



coefficients of its characteristic polynomial M.[x] 3 p{x; A) := 



Iax^ + IIa a; - HIa- 



Both for the expression of different models and for the construction of representation 
bases for redistribution, the following definitions are introduced 



-ij 



dxi dxi 

_ lA.. 
"2k 3 



6.... := 



2£ 



i(5 

3"y 



(4.3a) 
(4.36) 
(4.3c) 

(4.3rf) 



Sin '• — 



Pi. 
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I I dui duj 



dx, 
duj_ 
dxi 



dxe 



-i—rdui 



duj 

dxe 
duj 
dxi 



pu[u\Sji - puW^Sii 



Pi\j ■= - pKu'^^ji - pu'-u'ftu 



p(dov) 
ij 


P^j~ 




„(dcv) 
ij 


A, - 




p(dov) 


^=Ps., - 


2 p c 


p(dcv) 
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(4.3e) 

(4.3/) 

(4.33) 

(4.3/i) 

(4.3z) 
(4.3j) 



where (|4.3flP is the symmetric tensor of the second moments of the fluctuating ve- 
locities {—prij are the Reynolds-stresses) with trace Ir = 2k (|4.3al) and nondimensional 
deviatoric projection the Reynolds-stress-anisotropy tensor bij (I4.3cp . e\j^ (|4.3b|) is the 
symmetric dissipation-rate tensor (I3.1&P with trace I^c^) = 2s (j4.36p and nondimensional 
deviatoric projection the dissipation-rate-anisotropy tensor b^.. (|4.3cij) - Sij (|4.3eP is the 
symmetric deviatoric (in incompressible flow) mean strain-rate tensor, f2,y ( |4.3j| ) is the 
antisymmetric (hence deviatoric) mean rotation-rate tensor, Pij ( |4.3ffP is the symmetric 
rate-of-production tensor of the Reynolds-stresses (I3.ip with devi atoric projection pj'^^^^ 



(|4.3gD, the synimetr ic tensor Dij (|4.3fep is the complement of Pij ( Naot et al. 1970l 1973 ; 
Launder et aLl[l975 ) with respect to twice the symmetric tensor Pg.. (I4.3ip which is the 



production by mean strain with deviatoric projection Pj'^''^-' (l43il), and the symmet- 
ric tensor P^. . ( |4.3j| ) is the production by mean rotation which is deviatoricQ so that 
pLT^ :=Po./(E13- Obviously 



||Oi||l3i|[43j||43il 

*J Sij ^ij Si 



P — P 



2P, 



Q,ij 



Ip = Id = Ipg = 2Pk 
Ip. = 



(4.4a) 

(4.4&) 
(4.4c) 

(4.4d) 



By (j4.4cp Dij is better interpreted as the difference between production by mean strain 
and production by mean rotation. In (|4.31 14. 4p k (|4.3ap is the kinetic energy of tur bulence, 
£ (14.361) is its dissipation rate, and Pk (j4.4(jj) is its production rate. Following iLumlev 
(|l978l ) the invariants (14. 2p lib and Illb of the Reynolds-s tress-anisotropy tensor (|4.3cll 
play a major role in turbulence modelling. More recently ( Jakirlic fc Hanialial2002l ) the 
invariants Ilb^ and Illb^ were also considered. Both (|4.3ap and e\'^^ (|4.36|) are pos- 
itive semideflnite (their diagonal components are always nonnegative and they have 3 
nonnegative eigenvalues), so that they have the same realizability triangle, and we may 



t Pa.,, ^ -2pu'^u'fQ.u ^ +2pu\u'p.H = ^2pu!.u'^Q.,i, = +2pu\u'iP.ii -Pa.,, =^ Pa.,, = 
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define the corresponding Lumley's flatness parameters ( Lumlev 1978[ ) 

A-~l + 27111b + 9IIb 
:=1 + 27111b + 9IIb 



(4.5a) 
(4.5&) 



Both A (ji^^ and A^ ((436|) vary in [0, 1] with A — denoting 2-C turbulence (jSimonsen &: KrogstadI 
2005h . A — 1 denoting isotropic turbulence, Ag = denoting 2-C dissipation and Ag = 1 
denoting isotropic dissipation. 

4.1. Typical models for (pij 

It is beyond the scope of the present work to make an exhaustive review of existing mod- 
els. We concentrate instea d on models representative of wide ly used families of cl osures 
(lrr: 



Launder aLlll975l: G L: lGibson fc Launder Il978l: SSG: ISpeziale et aLiri99ll : says 



So et aLlll996t S: ISugal 120041: JEJK: Ijakirhc et a/.ll2007[) . These models are briefly sum- 



marised in the following, using the forms given in the original papersjj before proceeding 
to their representation in a common basis. 



4.1.1. iLaunder et a/.lll975l lrr 

The LRR model ( Launder et aLl[l975l ) is an early attempt to apply tensorial representa- 
tions for the closure of the order-4 tensor deflned by the volume integrals in ( 14JJ, leading 
to the same tensorial representation of 4>ij as the one which was obtained bv lNaot et al 



( 19701 . Il973h applying a closure for the 2-point correlations in (|4.ip . It reads 



'ij J LRR 

(s;2J), 

ij J LRR 

(r;W), 

ij J LRR 



,(s;u))i 



- 2cLRRi/oe6. 



( ^ r 4- 8 ^ p(dov) 

+ CLRji4/GL„ {-Pij -^ij) 



ij J 



(4.6a) 
(4.66) 

(4.6c) 
(4.6rf) 
(4.6e) 



Clrri — 1-5 ; Clrr2 — 0.4 ; Clrr3 — ; Clrr^ — ^ (4.6/) 

where Plf°^^ ( |4.3g[ ) is the deviatoric projection of the production tensor p. lap , d'"-^^^'^ 
(|4.3fep is the deviatoric projection of the difference between production by strain minus 
production by rotation (|4.4cp . The wall-echo terms (|4.6(jl I4.6ep do not use the tensor 
e-mGuj (cn being the unit normal-to-the-wall vector; cf i|4.1.2l) . but only a wall-proximity 



function /gl„ (4.7/), which is inversely proportional to the distance-from-the-wall n made 
nondimensional by the turbulence lengthscale ( |4.7/[ ), an d which discuss ed in further 
detail in M.1.21 Although superseded by more recent models ( IIaniali(j ll994') the lrr pro- 



posal is included here as a baseline reference, not only because it is one of the most widely 
known closures, but also because it can be used to highlight the problems encountered 
when using a linear model with constant coefficients. 

t Several models {eg S: ISugal |2004| : jejk: iJakirlic et al\ l2007f ) are written us- 
ing atj :— 2bij to denote an alternative anisotropy tensor (|4.3c[l . and the traces 
A2 := tra^ = aijaij = Abijbij = -8IIb and A3 := tra^ = aijajkCiki = Sbijbjkbki = 24111b. 
These substitutions were directly made in the model expressions. Lumley's flatness parameter 

(jLumlevlllQTi 'l is the same in both bases, viz A *W 1 -|- 27111b -|- 9IIb = 1 - §(^2 - A3). 
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4.1.2. iGibson fc LaundeilllQTSl . GL 



The Gibson-Launder model introd uces the modern form (IHanialic 
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for the wall- 



echo terms, by applying the ideas of Daly fc Harlow ( 1970l) and Shiij (|l973l ) to include 
both slow and rap i d wall-echo term s. It uses the isotropisation-of-production (ip) model 



( Naot et ali 



ig rap i Q waii-ecno term s, it uses tne isotropisation-oi-proauction [IP ) moae 
llQTOl : Revnoldj 1974 ) for the quasi-homogen eous part of r apid pressure- 



strain, along with the linear return-to-isotropy (ri) model of lRottal (|l951(a ) for the quasi- 
homogeneous part of slow pressure-strain, and reads 













._ „ p(dcv) 




' f ( 



,(s;to)i 



2 GL 



(4.7a) 
(4.7&) 
(4.7c) 

(4.7fl!) 



Vki jGLCnfcenj 

i^'-'^^l e e 



% lGLe„,e„, 



km 



(4.7e) 



Cgli = 1-8 ; Cgl2 = 0.6 ; c^^^i = 0-5 

3 

/gl„ — — ; as = 0.3 ; Kvk = 0.41 ; 1^ 



(4.7/) 



The tensor e„;e„^ (e„ being the unit normal-to-the-wall vector) was introduced by 
Shhl (Il973l). foll owing earlier attempts to modelling the surface i ntegrals in (|1.2I) by 
Daly fc Harlow ( 1970l ). to account f or the wall-suppressio n effect (jShir 1973 . p. 1329) 
which increases near- wall anisotropy (iDaly fc Harlow 1970l p. 2638), and has since been 
generalized in near- wall modelling (Hanialic 19941 ). The function /gl„ (|4.7ri| I4.7ep is 
a wall-proximity indic ator, which activates the wall-echo terms near solid boundaries. 
Daly fc Harlow! ( 1970l (20-22), p. 2638) have suggested that this wall-proximity function 
should be inversely proportional to the the distance-from-the-wall n made nondimen- 
sional by the turbulence lengths c ale £j , and this wa s also adopted in the defi nition of 
/gl„ ( |4.7y| ) bv lGibson fc Laundeil (|l978f ). As shown in lGibson fc Laundeij (|l978[ ) the con- 
stant in the definition ( |4.7y| ) of /gl„ can be related to the von Karman constant Kvk and to 
the value of the e quilibrium structural parar neter a,, := k~^ (— m'^ Q for simple shear flow 
( Townsendlll956l ). Althoug h the lrr model (Launder et al. 19751 ) includes an IP variant 
(lrr-ip), the GL proposal ( Gibson fc Laundeijll978 ) includes the general tensorial formu- 
lation for both slow (|4.7(ij) and rapid (|4.7e[) wall-echo ter ms, and has certainly been much 



more influential in modern near-wall modelling work ('Shimj 



19891: IShimalll993all6l: lHa,nialic fc JakirhcllT99 3: Hanialic 199^ 



1988: Launder fc Shim; 



Hanialic fc Jakirlic Il99 



Gerolymos fc Valletl[200l[|2002l:|jakirhc fc Haniaha2002l:IGerolymos et a/.ll2004l:|jakirhc al. 
20071) . 



4.1.3. iJakirhc et aLl 120071 JEJK 

The JEJK model ( Jakirlic et al. 2007^ is included as a typ ical representative of the 
family of quite successful models ( Shimal 1988t Launder fc Shinia.1989. : .Shima.l993a. ,fa 



28 G.A. Gerolymos, C. Lo, D. Senechal and I. Vallet 

Hanialic & Jakirlic 'l993; 'Hanialic 'l994'; Hanialic & Jakirlic 'l998!; Gerolymos fc Vailed 
200L ,2002: .Jakirlic fc Han ialic, 2002 : Gerolymos et aL.,2004) which have been developed 
using a quasilinear isotropisation-of-production (ip) tensorial formulation for the rapid 
volume part (/()|^'^\ with a corre sponding wall-echo term using a tensorial form 

based on a unit vector (|Shiij|l973l : iGibson fc Laundeiill978l ) pointing in the dominant di- 
rection of turbulence inhomogeneity. The dependence of the coefficient-functions of these 
models on the invariants of the anisotropy tensor (|4.2[ 14.51) is t he key element which 
greatly enhances th eir accuracy and ph ysical consistency ((Sambasivam et al\ 2004 ). Fur- 
therm ore, the JEJK ([Jakirlic et model, follows the suggestion of lHanialic fc Jakirlig 
(Il993l) to include not only the invari a.nts of ha (I4.3cl) but also those of &g (14.3 dll in the 



defin ition of the coefficient-functions ( Hanialic fc Jakirliclll993 . 1998 : Jakirlic fc Hanialic 
20021) . It reads 



[H'ij J JE.IK 

Pij JjEJK 
t>lj ''],IE,7K 



- pel 



, (s;i(;)i 



r ,(r;to)i _ (w) r 

[<Pij JjEJK ■— C]EJK2 /jE.IK 



(2c]E,iKi6y - 4c',E„^^(6ifc6fcj + |Hb%)^ 

„ p(dcv) 
CjEJKi /,IEJK„ ^^fe^rnSrjjj,,,;^ 6^,^,,,.^^ S- 



2 IV'/cj J,iE,iKe,,,j.,K^ 



(4.8a) 
(4.86) 
(4.8c) 

(4.8(i) 

7 

(4.8e) 



'-.1EJK2 



2.5A( min(^,-8Hb)) min 



"JEJKi ■ - 
KlEJKi 



V 



V 150/ 



k^ 



grad(A^T) 
|grad(A£,)| 



iin(/GL„,i) 



JEJKi 

3 



: — max | ttj yqCjejk, jg 



max ( 

1 CjEJKj 



64tt 1^ 
10 b' 2/ 



min(A|j) (4.8/) 



The JEJK (jjakirlic et a/.ll2007l) model uses a u nit vector e„,„,, (|4.8/|) p ointing in the dom- 
inant direction of turbulence inhomogeneity (jGerolvmos et a?.ll2004j) . and although the 
wall-proxir uity function /jrik... (|4.8/ i is based on the wall-proximity function /< 



(4.7il 



of the GL ([Gibson fc Laundeij 119781 ) model, and hence contains the distance- from-thc- 
wall n, it could easily be made fully wall-normal-free by appropriate redefinition of the 
wall-proximity function. 



4.1.4. ISpeziale et a/.lll991 



SSG 



The SSG (jSpeziale et a/.lll99lh model is based on general tensorial representation ap- 
proaches, and was developed as a closure for homogeneous turbulence. It is therefore 
interpreted here as a closure for the quasi-homogeneous volume part of redistribution 
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4.11) and reads 



mj \ssG 



(r;9J), 
ij J 



- CssGiPebij + CssG2 P^ibikbkj 

+ (cssG3 - 4sG3 \/-2IIb) Pk5y 

+ CssGiP^ibikSkj + Sikbkj — |lbs% 
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(4.9a) 
(4.96) 



(4.9c) 
(4.9d) 



Css 
C* 



3.4 
1.8 



4.2 



SSG3 



4 
5 

1.30 



CsSG4 
CsSGs 



1.25 
0.40 



(4.9e) 



The term ~c*g^^P]^bij (|4.9cp depends on mean velocity gradients because of the presence 

of Pir (I4.4d | 4.3<| and is therefore included in the rapid part 0-]^'^^ ( Shih fc Lumlev 1993; 
Ristorcelh et aLlll995l) . 



4.1.5. iSo et a/.lll99a SAYS 

As mentioned above ('M.l.'ll) the SSG mode l for i 



flows (jSpeziale et al\ Il991l ) 



was designed for quasi-homogeneous 
So et al\ ()1996t ) developed a near-wall version of the SSG 



model, by adding wall-echo terms to it. Incidentally, they used Pj;^'^^^ ( |4.3gP and D^'^'^^^ 



(j4.3fep for an equivalent with (|4.9cP tensorial representation of [4>'. 
The model reads 



(r;aJ)i 



[V'i j J SAYS 
J SAYS 
J SAYS 



,(s;9J)- 



SAYS + 



i (s;U!)i 



[0. 



aC'';^')! 



-CssGiPebij + CssG2 P^ibikbkj + |llb<5ij) 



,(r;»)i (|lli]|T33Elill * p , 

'ij JssG — CsgQ^-TkOij 

f- (cssG3 - |CsSG4 - CsSG3\/^2IIb^ pk5j. 

~ i(CsSG4 + CsSG5)Pij ^ 
4 l,LssG4 CssGs 



(s;u))i 
ij J 

(r;w) 



[Vii J SAYS JSJ 



1 14.1061 1 



A(^;a3) 



]ssG /si 



]sAYS :— — CgsQ^Pk^ij 



* p(dov) 
"says ij 
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80 



4(CsSG4 + CssGs) — ■ Q^SAYS 
4(CsSG4 ~ CSSG5) ='■ /3sAYSi — ^ 

a:,,, - -0.29 



CSSG3 



27.SAYS kS'ij j /sAYS„ 
2'ysAYSi 



jCssG, 



/sAYs„ := e y^oo) 
7|,,s = 0.065 
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(4.10a) 
(4.106) 



(4.10c) 
(4.10d) 
(4.10e) 

(4.10/) 



The wall-echo terms (|4.10d I4.10ep use a wall-proximity function /says„ ( |4.10y[ ) which is 
based on the turbulence Reynolds-number Rct ( |4.8J only. 
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4.1.6. ISugall2n04l Suga 

This closure for 4>ij is one of the recent representatives of the so-called 

two-component-limit (tcl) closures, which are based on the representation for the quasi- 



homog eneous part (^,-^^ developed by Fu (Il988l). to which i nhoniogeneous terms were 



added (|Craft fc Laundeij|l996l:ICraftlll998l:lBatten et a/.lll999l:ISugall2004HAl-Sharif et al 



these variants differing mainly on the particular choice for the free coefficient 
function Cg^ which multiplies in ()4.11c|) the tensor i^TCLi^ ( |4.11J and on the model for the 
inhomogeneous terms. A systematic study of the differences in results between variants 
of the TCL has not yet been conducted. The S (|Suga.2004) variant of the tcl mo d el use s 
a simplified representation for the inhomogeneous terms, following Batten et al. { 1999t ). 
It read^ 



,(s;tu)i 



s + 



,(s;to) 



]s + 



- 2(csi + c'sjpebij - 4cs^c's^pe{bikbkj + ^HhSij) 



0.6P.' 



(dcv) 



1.2Pkfe 



0.2- 



2S, 



ki 



(A) 
3 



--fsj\i^^didn{didj ~ ^dkdkSij 



(4.11a) 
(4.116) 



(4.11c) 



(4.11d) 
(4.11e) 



J^c.., :=(^ + 2IIb)if^^) 



^Pk(26,, -4(6,fe6fe, + |IIb%) 



10 i^'e^'i 



H '^k"-j 



10 



k2 



{Dek ~ Pik) 



(4.11/) 



d. 



(4.115) 



f In the present paper we conformed to the usual definitions of strain-rate (|4.3e|) Sij and ro- 
tation-rate ( |4.3y| ) flij tensor s, which are also adop ted in the more recent papers on the TCL (two- 
component-limit) model (eg'Al-Sharif et al''201(i). The equations arc identical wit h those of ear- 
lier p apers on TCL {Fu 1988; Craft & Launder 1996; Fu & Wang 1997;^ Craft 1998: iBatten et al\ 
[1999), including [Sugal IQ004), where [SijlrcL ■= 2Sij and [f7»j]TCL := Sfiij are used instead, with 
appropriate multiplicative modification of coefficients. 
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fA 



A 



A < -2- 

^ - 100 

5 ^ A ^ 1 
100 ^ ^ 10 

^ - To 



^si 10 



CS2 



™n<^ ^ l-cxp(- 



-A'^Re^ 
100 



32 ^ 
10^ 



= min(^, v^) 
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10 ■ 



4\/6 mil 



3IIIg 



3 + -(^y=8II^+y^8II^ 

Re-r - 55 ^ 



(-2IIg)i 



3(1 - \fA\ min (l,max(0,l - 



70 



lib (1 - VA) min ^1, max (0, 1 



i?eT - 50 ^ 
200 ' 



10 



3/s 



(4.11/i) 



The above equations are presented in their original form (|Sugall2004r ) and can be reduced 
by projection on an irreducible basis (c/ ii4.2p . Notice also that the presence of the 
invariants of S and J7 in the definition of Cg^ (j4.11/t[) . in particular the fraction which 
forms the second term of the three-term-sum expression of Cg^ (|4.11/tp . make the tcl 



model for [0i5''^'']s (|4.11cp weakly no nlinear with mean ve l ocity gradients, contrary to the 



^ \ 1 1 1 I — — 1 1 1 r 

usually applied linearity constraint (jLaunder et aLl 119751 : lLumlevlll978t iRistorcelli et 

I1995I) implied by 



4.2. Tensorial representations for including wall-echo 

To compare different models and analyse differences between them, it is necessary to 
project (rep resent) these models in a common tensorial basis ( Spencer fc Rivlin 19591 : 
Smithlll97ll) . This is obviously not the ca se for the models co nsid ered in the present 



as a linear combination of 



work ( ^4.1). Several of these m odels (lrr: Launder et al. 19751: GL: Gibson fc Launder 
1978t SSG: ISpeziale et "aLl ll99lt SAYS: ISo et aLlll996l) represent the closure for the quasi- 
homogeneous part of redistribution ^^J'^ (j>[j'^'^ + ^ij'*^^ 
the tensors bij (|4.3cp . Sij (|4.3ep . ( |4.3ffP , and Da (14.3/ip. and their p roducts which are 
linear in mean velocity gradients. The SSG model (jSpeziale et aZ.lll99l[) was written using 
the tensors bij (|4.3cp . Sij (|4.3ep . and ( |4.3j| ), and their products which are linear 
in mean velocity gradients, but can also be expressed using the previous tenso rs {bij, 
Sij, Pij, Dji), which are the ex pressions used in the SAYS model (|So et all 119961 ) . which 



uses SSG (jSpeziale et al\ Il99l[ ) for the quasi-ho mogeneous part of the model (I4.10cl) . 
Finally the models of the tcl f amily (Ei 1988; C raft fc Laundeij|l996l:lFu fc Wa ng 1997t 
CraftHl998t batten et al\\l99^ ISugall2004l ) use (|OT|) all these tensors {b,j, S^-j, fy , Aj, 



but also rij (14.3 al). and their produ cts which are linear in mean velocity gradients. 
It is established ( Ristorcelli et al. 19951) th at the most general model developed from 
( 4.1 all assuming homoge neity ( LumlevI 1978 ) is expressible in the irreducible polynomial 



Spencer fc Rivlinlll959t) representation basis of isotropic symmetric deviatoric tensor- 



valued functions of the 3 tensors bij (|4.3cp . Sij (|4.3ep . and VLij ( |4.3J , which is formed by 
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^tensors feo 'Ristorc elli et~ai. 19951 (42 j . p. 131), and t his is the classical represe ntation 
basis <B[b, S, U\ ( Gerolvmos et alMoiOm . Alternatively (|Gerolvmos et aLll20106h we can 



construct a basis of 8 symmetric deviatoric tensors Q3(p) [b, S, Pq; Pg], where we use 



explicitly Pq 



Vth, h^Vl - J7b^ and h^Vlh - bflb^, and Pg,, ^ -\^buS^ 



and its products with bij (|4.3cp to express bS + Sb 



fliibij^ and its products with bij (|4.3cp to express bfi 

2t _^ 2 



Siibfj — gIbs) ^ s^^ij 



3^bS^ 



Sb-" 



2 



1x2 a la- 



Passage matrices € 



from (to) <8[b,S,f2] to (from) <8(p) [b, S, P^; P^] are quite 



simple ( Gerolvmos et al. 201061 Tabs. I and II, p. 2). We will represent each model for 
the quasi-homogeneous volume-terms 



2 



(n)ij 



n=l 



{n)ij 



(4.12) 



where the slow part of the model e ventually (ym 
contains the dissipation anisotropy ( Lumlevlll97S 
tensors 



= 1) but not necessarily (xm^ = 0) 
and the symmetric deviatoric basis 



H 



(2)U 



H 



H 



H 



(5)»J 



H 



H 



(7)»i 



H 



-bigbgj + |llb 5ij 



-Pg — 



' pe ^ 



--bilPse, 



-bikbkiPnij + Pna^ekbkj 



(4.13a) 
(4.136) 

(4.13c) 
(4.13rf) 
(4.13e) 
(4.13/) 

(4.135) 
(4.13/i) 



were made nondimensional by scaling with powers of k (|4.3ap and e (j4.36[) . The represen- 
tation of different models for t he quasi-homogeneou s part of redistribution i n jB[b, S, fl] 
is av ailable in the literature (jSpeziale et all llQQlt iRistorcelli et al. \199A IFu fc Wand 



1997^. and it is straight forward to project these representations on *8(p) [b, S, Pfj; P5 



( Gerolvmos et al. 2010^, with the exception of the tensor F^chij (|4.11cp . which is invari- 
ably given in the reducible form ( |4.11J . Its expression in ^(p) [b, S, P^; Pg] is obtained 
after straightforwar d albeit leng thy calculations, which require the use of the C ayley- 



aiter straigntiorwar d. albeit leng tny calculations, whicn require tne use or the y ayley- 
Hamilton theorem (jRivlinl [19551 (4.24), p. 689) and its extensions (|Rivhnl [l955i (4.21) 



f These can be reduced to 7 using RE-projection (|Rivlin fc EricksenI 119551 . pp. 380-382) to 
express b^flb — bf2b^ as a li near combin ation of the other basis tensors, thus obtaining a 
functional representation basis (jSmithlll97lf ). ie a basis where isotropic tensor- valued functions 
are represented with coefficients which are rational (not necessarily continuous) functions of the 
invariants. 
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and (4.22), p. 688), as 



Ft, 



ps 



pe 



^Wv + —^H(2)ij + (f (^ + lib) - f 



|lIb-IIIb))i? 



+ (23 
15 



2iib)if(4te + (^ + 6iib)if(, 



5)y 



13 TT 



lIi/(7).,+4i/(8),, (4.14) 



Since in the present work we are interested in wall effects, we must also include inho- 
mogeneous terms in the representation. Curre n t and past modelling statu s (jShiiJ |1973 ; 
Launder et al. 19751 Gibson fc Launder 1978 : Gerolvmos fc Vallet 12001 ) usualljU in- 
cludes a unit vector pointing in the dominant direction where turbulence inhomogeneity 
is the highe st. In wall-topology-dependent models this unit vector is simply the normal 



to-the-wall (Shiijll973t Launder et 



topology-free models (jGerolYmos fc Valle' 



1975 



. .Gibson fc Launder 1 9 78| ) . e„. while in wall- 
20011 : iGerolvmos et al. US) this a unit vec- 



tor e*! determined from gra dients of the local state of turbulence (length and time scales 
and anisotropy invariants; IGerolvmos et al. 2004 (3), p. 410), or i n the context of e l 



liptic blending models by the gradient of the blending parameter dOceni et al. 2008 ). 



Several wall-normal-free models which contain i nhomogeneous terms (jCraft fc Launder 
19961: ICraftI [1998I : iBatten et aLlll999l : ISugal l 20o3 ) contain such inhomogeneity-direction- 
unit-vectors, although this i s not always stated explicitly, and in some ins tances 2 or 
3 different vectors are used ( Craft fc Launder 1996 : Craft 1998t Suga 2004). Let us for 
the time assume that the modelling strategy includes a single unit-vector. ISmithl (jl97ll) 
develops representations making explicit use of vectors, and applicati on of this approach 
leads to the functional representation basis in IGerolvmos et al. I (l2004l . (31), p.4 28), built 



from hij, Sij, J7y and Cj. Obviously this basis, written using explicitly Ci, can be rewritten 



using the deviatoric projection of the symmetric tensor Ci 



e*,, - il3 := 



I e„ — -lis geometric normals 



Ci <8) 6*1 - 



WNF models 



(4.15a) 



which, being the deviatoric projection of the tensor product of a unit vector (e,,) by itself 
satisfies 



2 14.151 , 1 r \/ 1 r \ ^ -* 




2t 



^ = 



n I 14.13at 1 9 9 I l4.13at , / , 



2T 

9^3 



I„ «^0 



llr 



iiii, — 27 



(4.156) 

(4.15c) 

(4.15d) 
(4.15e) 
(4.15/) 



As a consequence the basis should only contain elements linear in 77, since all higher pow- 
ers are expressible us ing (14 . 1 5 61 14 . 1 5 c| and relations derived from these. The functional 
representation basis (jCerolvmos et al ] l2004i (31), p. 428) is not very practical, because 
the repr esentation of, eg, the classi cal echo-term of the isotropisation-of -production rapid 
closure ( Gibson fc Launder 19781 ). cannot be represented explicitly ( Gerolvmos et al 



f As stated in ^ the alternative approach of applyinR ellip tic relaxation to the quasi-ho- 
mogeneous model to account for wall-blockage (|Durbinl 11993 !) has already been studied by 
[Manceau et al.. (|2001i ). and is not investigated in the present work. 
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2004"), th is being a rccurcnt problem with functional representation bases (jRistorcelli et al. 
1995.; Wallin &: Johansson 2000.,'). F or this reason, we use a polynomial basis, constructed 
following ISpencer fc Rivlinf ( 1959), who worked directly with tensors. The polynomial 
representation basis constructed ( Spencer fc Rivlin 19591 Theorem 1, p. 320) by adding 
■q (|4.15aP to the tensors generating Q3(p) [b, S, P^; Pg], viz bij (|4.3cp . Sij or Pg.. (|4.3ep . 
and Pfi.. ( |4.3j[ ), to construct Q3(p)[b, S, Pfj, 17; Pg]- Since the representation must be 
linear in mean velocity gradients, and because of (|4.156[ H.lScp . terms 



(a) nonlinear in S and/or $7 (ie terms containing [S]"^ [Tl]"'^ with rii + ^2 > 1) 
(6) nonlinear in rj 



are omitted. This leads to adding to the 8 tensors of the representation basis (|4.13p the 
tensors 
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■=Vij 




(4.16a) 


H(r,- 2)ij 






(4.166) 
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■=bubekVk3 + mebekbkj - f Ib^r, ^ij 




(4.16c) 








(4.16d) 
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TT 
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(4.16/) 




■■=mibtkP-s^, + PsJikVkj - fl^bPg 




(4.16<7) 


H(r,- S)tj 


■■=V^^PsJkJ + htPs,,Vk, - §1^^^* _ 


S^j 


(4.16/i) 


H{ri; 9)ij 


■.^mebekPfi,, +PnJekVkj ~ fl^bPf, 


Sij 


(4.16i) 
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H{r,-W)ij 


■■^miP^Jkj + baPn,,Vkj - 




(4.16j) 
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'■=r]iibe,kbkmP Sry^j + P Subikbkmrjmj — 




(4.16A:) 




■='niePstk^kmbmj + buhlrnP Smk''^kj — 




(4.160 
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■—ViebimbrnkPn^j + Pnnbimbrnkllkj " 


- . 5 - 

^ r,b-h 


(4.16m) 


* 


■—ViePnek^kmbrnj + bubimP n^^^Vkj " 


2T ^ ■ 
3\b^P, 


(4.16n) 



If more than one inho mogeneity-direction vectors are used [eg Craft fc Launder 19961 : 
Craftlll998l: JSuga"2004|)_ then one ca n either repeat (|4.16p for each 77, or attempt a RE- 



projection (|Rivlin & Ericksenlll955l pp. 380-382) of aU 77's on a selected one, using as 
a second tensor, eg b, but also keeping in mind that the projection may fail if the 
2 selected tensors are not linearly independent, for the particular case of interest. In 
practice only a few of the tensors in (I4.16P are used {viz H(^j^. i^j^, 2)ij, 4)ij, 
■^iv; 5}ij ^^'^ P^im models for wall-echo considered in the present work ( §4.ip 
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Table 3. Coefficients of the tensorial representation of redistribution closures (|4. 131 14715] 14.161 14. 17p (for [<^ij '^^]M l|4.12|l only coefficients which are 
7^ for at least one of the mode ls are tabulate d) , for the mo dels co nsidered in the present work (lrr: [Launder et oZ.lll975l : GL: lGibson fc LaunderlflQTSl : 
SSGi lSpeziale et fflZ]|l99ll : SAYS: ISo et aZ.lll996l : S: ISueal [20041 : jejk: I Jakirlic et a0 [2007) with reference to the coeffici ents of the ori ginal representations 
( i]4.1|l . and indicator Xm^ (|4.12|I of the inclusion (xu^ = 1) of the anisotropy of dissipation in the slow volume part l|Lumlevlll978f ) of the model (recall 
that the SSG model is quasi-homogeneous and that the says model is its near- wall extension). 
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can be represented as 



n=l 



n=l 



As;w) 
ij 

14 



(4.17) 



where a different tensor i? 



H 



(1?;1 



, based on the unit vector (j4.15ap , was used i n the 



rapid term, to accommodate the different inhomogeneity direction used in the S (jSuga 
2004[ ) model (|4.11ep for the slow and rapid terms. 

All the models studied in the present work f H4.ip can be represented by l|4.12[ BTT7)) . 
with coefficients (Tab. [3]) which are readily obtained by straightforward co mputation from 
the ori ginal form of the model. Alternative equivalent representation bases (jGerolvmos et al 
201061 ) can be used, th e coefficients in these di fferent bases being easily computed using 
basis-change matrices ( Gerolymos et al. 201061 ). 



4.3. A priori assessment of typical models for (pij 

A priori assessment of the various models for the rapid pa rt of pressu re-str ain redistribu- 
tion (Fig. [T5|) indicates that the two more recent models (s: Sugall200"4 : JEJK: Jakirlic et all 



2007^ give the best estimate of ^ 



both globally and for the individual volume 

and wall-echo (^Ij'™^) terms (Fig.[T31). The prediction of the (f>yy component is the most 
unsatisfactory, both these models (s and JEJk) predicting the wrong sign ( Fig. [I^, in- 
dicating that the much more complicated tcl representation used for the S ( Sugal 2004 ) 
model is not sufhcient to correct this discrepancy, common to all of the models studied 
(SjI, with exception of the GL ( Gibson &: Laundeilll978[) proposal. Nonetheless, even 

in this case, the separate assessment of the volume f^^d wall-echo ((^[j''"'') terms 

(Fig. [T51) shows unsatisfa ctory prediction of the yy-component for all models, the GL 
( Gibson fc Laundeil Il978l ) closure overestimating wall-echo by an order of magnitude, 
and thus overcorrecting t he equally ov erestimate d volume term (Fig. | 13[) . In general, 
both legacy models (lrr: iLaunder et a/.i,1975, : GL: iGibson fc Launder 1978|), which have 
constant model coefficients (not functions of the invariants), often exhibit an order-of- 
magnitude difference with respect to the DNS data, especially when examining individu- 
'■^'^^) and wall-echo ((/'■^'"'') terms (Fig. [T31). Nonetheless, observation 



(Fig. [T^ shows a fund amental difference be- 



ally the volume 

of the shear component of wall-echo 

tween the wall-echo approaches, of GL ( Gibson fc Launder 1978[ ) based on the general 
tensorial formulation ()4.7e|) , and of lrr ( Launder et all 1975h which simply uses a wall- 
proximity function approach (I4.6e|) . The directionality information present in the GL 

closure predicts the correct form of (l)xv" \ contrary to the LRR 
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Launder et aLlll975l ) model (Fig. [T31) . The quasi- homogeneous SSG (jSpeziale et aLlll99ll) 
largely overpredicts and the corresponding near-wall model ('SAYSt lSo et al. 19961) 

has to introduce overpredicted wall-echo terms to correct this discrepancy (Fig. [T3)) . 

When assessing the models for the slow part (Fig. [T4)) we included in the DNS data 
the pressure-strain redistribution tensor p.2p augmented by the anisot ropy of dissipa- 



tion (|4.4(jj) . 4>ij — 2pebe.-, since several closures follow the suggestion of Lumlev ( 1978() 
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Figure 13. A priori com parison of various model s (lrr: 
iGibson fc Launde/ ll978l : SSG: ISpeziale et al\ llQQll : SAYS: ISo et al\ 
Dakirlic et a/,.2007; ) for the rapid pressure-strain redistribution 



Launder et al\ 119751 : GL: 
19961 : S: ISugal |2004 JEJK: 
0iJ'^'+</'ii with present 

DNS results for the volume (/"ij'^' and the wall-echo ^^J'™' terms, in wall units, plotted against 
the nondimensional distance from the wall (in fully developed incompressible plane channel 

flow (j)yz — 4>zx)- 



, (r;m) 



to include the anisotropy of dissipati on into the m odel for the slow part of 4>ij- This 
suggestion was motivated by the fact (jLumlevI 119781 ) that, in the context of SMCs, both 



and hp 



are modelled by a quas ilinear combination o f 6y and bjmbm-i ■ 



11 — - -cij — "J — 1 — — •-. — -tj ... i7n-"'[ ■ jTl\^6jj (|4.121 

Kl^ . Ob viously fFig. UM the LRR ([Launder et a/.lll975[ ). GL (jGibson fc Laundeijll978l ). 
and SSG ( Speziale eit aZ[ 199lh proposals include the anisotropy of dissipation into the 



model for the slow volume part 



2pebs^., while the JEJK (jJakirhc et al\ 120071) 



closure models 



alone. Examination of the S (jSugal 12004) closure shows that in 
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s (lrr: 


Launder et al\ 1 19751 


ISo et al\ 


19961: S: ISuEal |2004 



iGibson fc Launde 

iJakirlic et aZ.ll20d7l ) for the slow pressure-strain redistribution (f>^^^ = 9!>,-°'^' + (pij'™^ 
DNS results for the volume ^i^''^' and the wall-echo ^i'^'™'' terms, and slow redistribution aug- 
mented by the anisotropy of dissipation <j>l^^ — pe'^''^'^^', in wall units, plotted against the 
nondimensional distance from the wall (in fully developed incompressible plane channel 

flow (j}yz — 4>zx)- 



the buffer and outer regions (y+ ^ 20; Fig. [M]) it overpredicts (pf^ but closely fits 

— 2pebeij, while near the wall [y^ ^ 10; Fig. [Ti)) it closely fits (j^l^^^ alone, as it 
respects the two-component-limit (tcl) for A ^ Q. This behaviour is particularly visible 
for the xs-component (Fig. [14]). Concerning the global prediction of (jif^ (Fig. [14]) the 
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s (lrr: 


Launder et all 1 19751 


ISo et al\ 


19961: S: iSuEal |2004 



GL: 
JEJK: 



Jakirlic et 0^.11200% for the pressure-strain redistribution cf>ij with present DNS results for redis- 



tribution 4>ij and redistribution augmented by the anisotropy of dissipation (pij — pe\'^''^''^^ , in 
wall units, plotted against the nondimensional distance from the wall j/"*" (in fully developed 
incompressible plane channel flow <j)yz — 4>zx)- 



JEJK (iJak irlic et "all l2007l) closure gives the best results, followed b y SAY S feo et a/J ll99l 
3 ISuga ^ 



and S (jSug a 20041. Notice that the global prediction of the SAYS (jSo et a l. 199Q) closure 

(s) 

fits alone, while its volume part which is given by the quasi-homogeneous SSG closure 
fits (/>j-!'^^ — 2p£beij ■ a consequence the wall-echo term 



as cancelling out the dissipation anisotropy included in 
(jSpeziale et aLlll991l ) closure, and to a lesser extend the lrr ([Launder et 



ffij J say; 
,(s;aj), 

Jij I SSG- 



must be interpreted 

Notice th at the ssG 
T975I) 



and 



GL ( Gibson fc Laundeilll978[) models, perform reasonably well in predicting 0-^^ — 2yoe6e.^ , 

and not (t)f^ (Fig. [HI). 

Concerning the global prediction of (Fig. [15]) the trends observed in the detailed 

predictions of 4? = ^J'""^ +i£!L(Fi£JIl and (/.g^ = ^g'*"^ + (^^8- ^ a^*^ 



confirmed, in that the JEJK (Jakirlic et al. 



by the S ([Sugall2004[ ) and SAYS ([So et al 



■20071) closure gives the be st prediction, f ollowed 
19961) models. The SAYS ([So et a/.|[ 19961) mode l 



does a very good job in correcting the quasi-homogeneous SSG ([Speziale et al\ 119911) 



closure (which also contains the anisotropy of dissipation) for near-wall turbulence, but 
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closu res whose vol ume part does not include the dissipation anisotropy 
JEJK: 



'2pebeij (such as 



Jakirlic et al. 2007. . and to a lesser extent S: lSuga,2004) are much easier to adapt 



to near- wall turbulence modelling (Figs. [TMTSt . 

Notice, however, that even the best models, which predict satisfactorily the behaviour 

(tog 



of (j)xx (in plane channel flow u'u' is produced by Pxx — —2u'v'dyu), are less accurate in 
the approximation of the redistribution for the two other normal components, especially 
in the near-wall region (?/+ ^ 10; Fig. [T4|). 



7(P) 



4.4. Typical models for d^^ 

There are re latively few models for pressure diffusion which were r ecently reviewe d in 
Vailed (l20Q7|). We examine here most of these m odels (|Hirtl (Il969l) : iLumlevI (|l978t ): |^ 
(|1993[ ): DRDL: lDemuren et aLl[l99l I Vailed (|2007t )). Notice thatlSuga| (|2004l) has proposed 



a model for the rapid part of pressure diffusion d^^''^\ indeed for the rapid wall-echo 
but the proposed model is identically in f ully developed pla ne channel flow. 



and cannot therefore be assessed in the pres ent work . IStraatman (Il999l) has studied the 
calibra tion of the coef ficients of the model of iLumlev (1978) ( ^4.4.21) . using the tensorial 
form of lLumlev (1978) to model the triple velocity correlations, and hence belongs to the 



same family as lLumlevI ()l978l) and I Vailed (l2007l) (for d 



)• 



4.4.1. lHirdll969l 

The early attempt of Hirt ( 19691 ) at modelling the slow volume part of dm El 



is a tensorially in v ariant gene ralization of the specific to boundary-layer flow proposal of 
DonaldsonI (fl969l) . iHht l ("1969^ has no model co efficient (equivalently the model-coefficient 
is taken to be 1), while Daly fc Harlow (Il970i) who used the same formulation include 
a model-coefficient, although both [Donaldson (|l969t ) and Dalv fc Harlow ( 1970f ) set the 
corresponding model-coefficieniU to 0. This model reads 



[O-ij J Hirt 



rj(p;s;aJ)i 

["ij J Hirt 



r ,(p;s:«;)i 
[Uij J Hirt 

r^(p;r;2J)i 

J Hirt 

r ,(p;r:iu)i 
I'^ij J Hirt 




\ Hirt 



CnirtP 



,{p;s;w) 



(p;naj)i 




(4.18a) 

(4.186) 

(4.18c) 
(4.18d) 
(4.18e) 



Cnirt — 2cLunilcyCDH ] Cl 



umlcy 



22 
100 



Had we used, a model-coefficient Cnirt 

equal to 1 instead of ( |4 18jD 



(4.18/) 
as used by iHird 



(|l969h and suggested by I Dalv fc Harlowl (|l97fl[ ). the values for d,^^'' would have been an 



order-of-magn itude too high compared to DNS data ()Kim et aLlll987H Moscr et al. 199 



Hovas fc Jimenez 20061: iGerolvmos et aLll2010cri . To compare the early (Donaldson. 4 9691 : 
Hird 1969t Dalv fc Harlowl 1970 ) tensorial form (|4.186p with later pro posals we used the 
same coefficient (4.18j) as the one (4.21^) used in the model of |pu| (1993), which has 



a similar although not identical tensorial representation (|4.216|) . but with the opposite 
sign. 



t (|DonaldsoDlll969l . (11,12), p. 273) and Sdh (|Dalv fc Harlowlll970l . (16), p. 2637) 
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Lumlevl (|1978l ) made analytical calculation in Fourier space to obtain a quasi-homogene- 



ous expression for the slow volume part of the pressure- velocity correlation p'u'^, whose 
validity is limited to weakly inhomogeneous flows 



1 



[P "iJLumley 



(4.19) 



ij jLumlcy -^["ij jLumlcy 



Lumlcy 



Lumlcy 



(p;s;2J)i 



Q I 

Immlcy : = (^CLumlcyp(MiUm'"m'^j^ + 



Lumlcy 

(4.20a) 
(4.206) 



r ,(p;s;tu)i 

["ij jLumlcy 


:= 


(4.20c) 


rj(p;r;ffl)i 

nj J Lumlcy 


:= 


(4.20d) 


r ,(p;r;™)i 

["■ij jLumlcy 


:= 


(4.20e) 



1 

CLumley — 5 



(4.20/) 



In the present a priori assessment the proposal (|4.19p of iLumlevi (jl97a ) was evaluated 
using the DNS-co mputed triple velocity correlations, and not a model for u^Ujui, a.s in 
previous stu dies (|Fulll993t IStraatmanlfl999t iDemuren eraZlll996t ISauret fc V^Ieti 120071: 
Valletil2007l ). Hence the model tested is not a SMC, but a third- moment closure (tmc) 
instead. 



4.4.3. lFul ll993l 

Ful (I1993I) used the proposal (|4.19p of lLumlevI (|l978f ) combined with the lDalv fc Harlow 



(| 19701) model for the triple velocity correlations 
of the pressure- velocity correlation p'u'- 



u'jU'jU'j^, to model the slow volume part 





:=[4-^^]lu + [4r^'"^]pu + Mif ^^^^]pu + id^-'^-x. 


(4.21a) 




d I k / 9k , , 9k \ \ 


(4.216) 


^d^'^'^Xn 


:=0 


(4.21c) 


["■ij Jfu 


:=0 


(4.21d) 




:=0 


(4.21e) 




Cpu — Cnirt — ^CLumlcyCDH 


(4.21/) 



The model has been applied to the study of the plane/round jet anomaly (Pu Il993l) 



4.4.4. iDemuren et a?.lll996l DRDL 

Contrary to the previous proposals, IDemuren et al. (1996) include a model for the 
rapid-part of the pressure diffusion term. The slow part was modeled by using a non- 
local elliptic relaxation approach applied to a baseline local (single-point closur e) model 
for the slow pressure transport p'u[, obtained by applying the model (|4.19p of iLumlev 
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( 197ah to the lMellor fc Herrind (|l973h closure for the triple velocity correlations 



k2 /9« 



dxk 



dx,. 



dxi 



^Lumlcy [^i^m'^^m] MH 



l4:T8/t 1 

CLuinley — 5 



Cmh. 



i Chl. 



Chl. 



11 
100 



(4.22a) 

(4.226) 
(4.22c) 



The lMellor &: Herrind (Il973f) clo sure for the triple velocit y correlations (I4.22a l can b e ob- 
tained ( Demuren fc Sarkaill 19931 ) by isotropisation of the Hanialic fc Launder ( 1972 ) clo- 
sure, with the relation (|4.22cP between the corresponding coefScients C mh. and Cht,.- The 
initia l estimate of pressure diffusion using (j4.226p in a plane mi xing layer ( Rogers fc Moser 
1994h moti vated the non-l o cal el liptic relaxation approach ( Demuren et all 1996f l. de- 
scribed bv IPemuren et al. ( 1994 ). The closure for the rapid part of pressure-diffusion 



models directly d^^''^''^-' as being proportional to the mean velocity gradients. 



[d. 



ij I 



■.=[d. 



(p;s) 



(p;s;<»)n 



[d. 



(p;r;<»). 



(p;r;ii?)i 



(4.22d) 
(4.22e) 



(p;s;u>)i 
ij J 

{p;r;aj)n 
ij J 



d 

^ I [P'u'illMHdje + [p'Uj 
(p;s;2J) 



(p;s) 



[p'u'i]i%TJSje + [p'u'jli^fiiUie 



-dUr. 



:=0 



(4.22/) 

(4.22<7) 
(4.22/1) 
(4.22i) 



wher e ERd is the elliptic relaxation operator of lDurbinI (|1993l ). defined as ([Demuren et al\ 
19941) 



L^V^ER„(/)-ER„(/) ; ER„(/) 







Lo := CDi^Tmax(l,CD„i?e. 



r - -1 



Co =70 



(4.22j) 
(4.22A;) 



Cdrdl,, 



J_ 

10 



(4.220 



In the a priori tests the operator ER^ ( |4.22ffP was computed by solving nurn e ricall 



the elliptic equation (|4.22ffP using 0{Ay^) space discretization. iDemuren et ali (|l99' 



recommend Cdrdl^ G [jqj yo]i and we have chosen the lowest value 
because the model generally overestimates d[^''^''^^ (of 



m 



To (ESSI 



4.4.5. Vallet 2007 



Valletl (12007|) suggested a model for pressure-transport p'u'^ based on an associated 



model for the triple velocity correlations which contains both slow and rapid terms, 
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following the analysis of lYounis et al. ( 2000l ) 



— \U.U -U.: 



- — \ u- u ■ u; 



allot,, ^ (^^1™,„^J^ + S^jmm'^*^ + ^'Imm'^u) (4.23a) 
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where the slow terms are based on the tensorial form of iLumlev ([l92i) 



du'iU'f 



dxi 



(4.236) 



with the coefficients Clsb^ and Cls b„, (|4-23./|) op timised by ISchwarz fc Bradshaw ([1993)0 
and the rapid terms introduced in IValletl (I20Q7D are based on an analogous tensorial form 



(r) - du'^u'g 



dxk 



Jkj— H Jkr 



dxk 



dxk 



(4.23c) 



By analogy to the relation (|4.19p of iLumlevI (|1978l ). slow and rapid pressure-transport 
terms are modelled as 



p "iJvallet 
p [P "iJvallct 



\ ' t / l(^) 
Cvallctsp Pi^m'^^mJLSB 



Cvallotsp Wi^'m'^'n^vlllet 



(4.23d) 
(4.23e) 



The model coefficients are (jValletj 120071 ) 

Clsb,, =0.0980 ; Clsb,, =0.01265 

Cvallct,, =0.0001 ; Cvallct,, = 0.0001 

cvaiiew :=0.085 (l + min (o.5, A^^^^-^^^ni-^^)) 



(4.23/) 



Notice that Cyaiictsr. is functio n of the flatne ss parameter A (|4.5aP to accommodate the 
model to wall turbulence. The IValletl (2007) model reads 



i-^lvallct 



_rj(p;s:2J)i 



(p;s;2J) 
ij 

(p;s;ui)i 



1-3 
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(p;s;iu)i 
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(p:r;9J)i 
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Jvallct 



Jvallct 



d 

dxi 



(s;QJ) , 



J vallct 



:=0 



(p;r;aj)i 



Jvallct 



l(r;W), 



[d 



(p;r;tu)i 
ij 



Jvallct 



(4.235) 
(4.23/i) 
(4.23i) 
(4.23j) 
(4.23A;) 



t lsb: ILumlevI ijigTST ) and ISchwarz fc Bradshawl l| 19941 ) 
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Fig ure 16. A prio ri comparison of vari ous mod els (H i rt: iHird 1 19691 : Lumley: iLumlevI Il978l : 
Fu: IFuI 119931 : drdl: iDemuren erallll996l : Vallet: IValletl 12007^ for the slow pressure diffusion 



^(p,s,w) ^gj.jjjg^ ^oii units, plotted against the nondimensional distance from the wall y"*" (in 



d)-' ' + d)-' ' with present DNS results for the volume d;?' ' and the wall-echo 

*j ''J ^ '■J 



fully developed incompressible plane channel flow d^' 



: d 



d' 



d> 



4.5. A priori assessment of typical models for d. 



(p) 



With the exception of the proposal of Hirt ( 1969I) . all of the models for the slow part d] 



7(p;s) 



are based on the proposal of iLumlevi (|1978[ ) that pressure-transport p'u[ is proportional 
J, ie to the transport of fluctuating turbu lent kine t ic ene rgy. The tmcQ a 



priori assessment (Fig. [16]) of the proposal (|4.19p of iLumlevI (|1978l ). shows not only 
that the coefficient CLumicy 



I ( |4.20j[ ) is too high with reference to the DNS data for 



^(p,s,Q3)^ but also that although the sign is correct for the shear component d^y^''^\ 

it is incorrect for dyy'*''^'' (Fig. [TB)) . A modification of th e coefficient Cmmioy ( 4.20^ in 
(|4.19p . either using a smaller constant Straatman (Il999l ) or a function ( |4.23j[ ) of the 
flatness parameter A ( Vallet! 2007 ). can improve the predicted magnitude_of d^^'*''^'' 



{eg the IVallct 200^ model in Fig. [T6]), but not the problem of sign. Hirt ( 19691 ) uses a 



different tensorial formulation (I4.18&[) resulting in the opp osite sign for compared 



to what is usually obtained from the relation (|4.19p of iLumlevI (|l978|j . but this does 



not solve the problem since for the model of Hirt ( 19691 ) the sign of dlfy^''^^ is correct 



but that of d^xy^''^^ is not (Fig. [16]). The SMC closure of P^^.t^-^u'^ using a relation of 
the form (|4.19l) . involve s a m odel for t he triple velocity corr elations, which influences 
the final prediction. 'Fu' fl99^ uses the jPalv fc Harlowl (|l970[) model for triple velocity 
correlations (|4.21&1 |4.21j| ) combined with the standard value of CLumicy = ^ ( 14.20Jj ), 
and the underestimation of the y-gradients of the triple velocity correlations by the 



J third-moment-closure using the DNS-computed triple velocity correlations u'jU'Yk in (|4.19[l 
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Daly fc H arlow mode l yields a corre ct final level for (Fig.fTH]). On the other 



hand, Vallet (2007.) uses the Lumlev ( 19781 ) tensorial formulation for the slow part of the 
triple velocity correlations (I4.23ap with the optimised coefficients of lSchwarz fc Bradshaw 
( 1994[ ). which yields slightly higher gradients for the triple velocity correlations, and 

(Fig. [T6)) by using a lower variable coefficien t (|4.23J 



recovers the correct level of d. 



(p;s;aj) 



function of the flatness parameter A in (|4.23ril) . Finally the model of iMellor fc HcrrinsJ 
( 19731) for the triple velocity correlations combined with the standard model (|4.19l) of 



Lumlev ( 1978t ). used for the slow volume part (|4.22a| - |^22cl) of the drdl ( Demuren et al 
1996[) model [d^^'^''^'']oRDL ( |4.22y| ), yields a very good prediction for the shear component 
d^xy^'^"^ but a very inaccurate one for the normal component dy^'y^'^^ fFig. [T6|. The wall- 



echo term \d, 



(p;s;u))i 



(|4^22ffD added by the elliptic relaxation approach (|4]22^^422| in 



the DRDL (jPemuren et a/.lll996i) model is designed to damp gradients. As a consequence. 



it adds a strong wall-echo term dlfy^'^^ ^(p.^.^j) 



to dyy ' improving the global prediction of d. 



(p;s) 
yy 



(Fig. [T5|. but adds an equally strong wall-echo term d^x-if'"''^ to the initially satisfactory 



volume term dxif' ' thus deteriora ting the global predic tion of dxy^' (Fig.fTB)) . 

Only two of the models (drdl: iDemuren et al\ Il996t Vallet: IValletl l2007t ) include a 
closure for the rapid volume term d- j'''^'^'' (Fig.[T7]). The levels of dl^'"^^ a priori predicted 
by the drdl ( Demuren et aLlll996f) mode l (l4.22/i| are an order-of- magnitude too high 
(Fig. [T7)) . The a priori assessment of the IValletl i2007h model 1^?^ gives reasonable 



(p;r) 



but the predicted normal component d^y^^ has the wrong sign 



magnitudes for d^^ 
(Fig. [13. 

A priori assessment of the global prediction of p ressure diffus ion d ffl {^}^- ^y 
the various models (g^l) shows that the models of l Valletl ()20Q7f ) and'0 (|l993l) give 
the best prediction. The drdl ( Demuren et al. 19961 ) model is incorrectly dominated 
by its rapid part (Figs. [16HT8| d,-j''^\ which is responsible for the global unsatisfactory 

prediction (Fig. [TBI), "^h e a priori assessment of the shear component d^xy shows that the 
models of Valleti i 2007 ) and Fu ( 1993 ) are quite satisfactory, except very near the wall 
(2/+ ^ 2; Fig. [TSj). Notice nonetheless, that all of the models for cjixy (Fig. [TS]) make the 
same error with the opposite sign in the near-wall region, so that the global prediction of 

(v) 

-\- dxy is not affected by this discrepancy. On the contrary, the prediction of 



'-'xy 



mi 

dm (Fig. [HI) is unsatisfactory, and although the models of lValletl (|2007[ ) and|^ (|l993l) 
are those having the smallest error, they yield, like all other models, the wrong sign for 
the normal component d'^y (Fig. [T5|) . 



5. Conclusions 

In the first part of the present paper (^J2I and fJSj) we studied the influence of wall- 
blockage on turbulent pressure fluctuations p', in fully developed incompressible turbulent 
plane channel flow, using DNS fi j2.ip . 

The instantaneous pressure fluctuations p' were separated ( ij2.2.1l and !j2.2.2l) into 
rapid, slow, ans Stokes terms (p' = p^^.^ "'"^'(s) ~^ p'{t)) solving the corresponding 1-D 
modifled Helmholz equation (j2.6p for each parallel- to-th e -wall wavenumber, using the 
standard Green's function approach (c/ Appendix) of lKim ( 19891 ). The new results in the 
present work concern the identification of the infiuence of wall-blockage (i j2.2.3l) . To this 
end we define the quasi-homogeneous volume part of the slow and rapid fields (p'^j-.^j) and 
P'{s-<a) ) solution of the corresponding 1-D modified Helmholz equation (|2.6ap , with 
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the same source-terms but with different, freespace, boundary-conditions (boundedness 
at oo), using the appropriate freespace Green's function (Tab. [2]). The wall-echo fields are 
then identified as the difference of the complete and volume fields {p'(^j..^'j ■= p' — Pji- qj) 
and 

The theoretical analysis (c/ Appendix) of the Green's functions for the generic 1-D 
modified Helmholz equation (|A 1[) was applied (i i2.2.4p to the determination of the error of 
the approximate method-of-images which replaces the exact boundary-conditions (|2.66p 
by using mirror images of the channel sources, one for each wall, demonstrating that the 
method-of-images is a high-wavenumber approximation. As long as the energy-containing 
wavenumbers of the source-terms are sufficiently high (equivalently the parallel-to-the- 
wall structure sizes f„ ^ 46, 5 being the channel halfheight) the method-of-images ap- 
proximation is satisfactory, and in this case f ii2.2.5p the interaction of wall-echo between 
upper and lower walls is negligible. When the relevant structure sizes are within this 
range (£« ^ 4(5), the wall-echo term, at the wall, is approximately equal (and of the 
same sign) as the corresponding volume term. These results are obtained by studying 
the Green's functions, and as such are independent of the particular distribution of the 
source-terms, depending only on the wavenumber. 

This methodology for p'-splitting {p' = p'^^.^j^ + P[r-tu) + P[s-<a) + P[s-w) + p[t)) ^as 
applied to DNS computations of low- Reynolds-number {Rct^ ~ 180) fully developed plane 
channel flow (Tab. [Ij. The results for the rms- values of the 5 p'-fields verify, near the 
wall, the approximate equality of wall-echo and corresponding volume terms ([P(r.rrj)]rms ^ 
[P(r;tu)]i-ms and [^(s.sjj) ] rms ~ [P(s;uj) 1 ims for y+ ^ 3). The rms-values of the wall-echo fields 
([^'(r-iu)]"ns and [P(g.„-,]rms) dccay with increasing distance from the wall, always remaining 
much higher than the Stokes pressure ([p(^)]rms)- The contribution of the 5 p'-fields to 
the correlations containing the fiuctuating pressure, viz pressure diffusion (E31), 

pressure-strain redistribution (pij (i i3.4|) and velocity /pressure-gradient 11^^ — <pij + d[^^ 
( ^3.5p shows that the approximate equality of wall-echo and corresponding volume terms, 
near the wall {y^ ^ 3), holds also for d\^\ (pij and 11^^13 

Examination of the relative importance of each of the 5 p'-fields in the pressure-strain 
redistribution tensor ipij m3A\i indicates, that for the (p^x component the rapid volume 
term is the dominant contribution near the wall (y^ ^ 10), while the slow volume 

part (/>j-j'^'' is the dominant contribution in the buffer and outer regions (y+ ^ 20). The 
opposite trend is observed for the (f>zz component, (plj being the dominant contribution 
near the wall (y^ ^ 10) and '/'Ij'^' being the dominant contribution in the buffer and 
outer regions (y+ ^ 20). Finally, the Stokes term (t>\^\ which is generally small, is the 
dominant contribution to the shear component (l)xy, near the wall ^ 8), indicating 
that p'^^^ is well-correlated with wall-shear. These results demonstrate that the relative 
importance of the rms-values of the 5 p'-fields does not determine the relative levels of 
the contribution of each field to 4f (gSH), (^33) and U^j f ^XS]) . 

In the second part of the paper (iJH) we use the DNS database generated in S|3]for the 
a priori assessment of representative models for 4>ij and d[^"^ . 

The models studied for the pressure-strain redistribution (SHI) w ere chosen to be rep 



resentative of modelling practices from the early developments (lrr: iLaunder et aLlll975 



f Notice that results using the approximate method-of-images should not be considered as 
a proof of this approximate equality, which is in-built in the method-of-images approximation 

(223}. 
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Gibson &: Laundeij llQTSi ) through the more gen eral q uasi-homo geneous tensorial rep- 



resentation appr oaches (SSG: Soeziale et aLlll99ll : SAYS: So et al. |1996.) to more recent 
pr oposals (jej k: Jakirlic et al. 20071 ). including a recent representative of the tcl family 



(S: ISugall2004l) . To provide a common representation of all these models we proposed 



( ^4.2p a tensorial representation including a tensor ij (|4.15ap . built as the deviator of 
the tensor product of the unit vector (pointing in the dominant direction of turbu- 
lence inhomogeneity) by itself, a long with the 3 generating tensors (b, S Cl) used for the 
quasi- homogeneous case (iRistorcelli et aLlll995t ) . The polynomial representation basis de- 
veloped in the present work (! j4.2p adds 14 new basis-tensors (|4.16p to the 8 tensors (|4.13p 
building the classical quasi-homogeneous basis. Furthermore, the rapid terms (containing 
mean velocity gradients) can be interpreted in terms of production by strain Pg.^ (I4.3zp 
and production by rotation P^.. ( |4.3j] ) , whose difference is the tensor Dij := Pg,^ — Pg 



J4.3/t[ 14.3 cP appeariiig in t he original formulation of several models (lrr: iLaunder et al 
[1975; SAYS: ISo adll99l S: [Suga 2004). 

The a priori asses sment of the different models for i^a (! j4.3D indicates that the jejk 



(Ijakirlic et al. 2007 ) model gives the best overall performance, followed by the S ( Suga 



200J) model. The design p hilosophy of these two models is fundamentally different. The 
JEJK ( Jakirlic et aLl 120071 ) model (Tab. [3]) uses few of the representation basis tensors 



(based on the ideas of isotropisation-of-production and of both slow and rapid wall- 
echo terms with inhomogeneity-directionality information) with coefficients which are 
carefully optimised functions of the invariants (both of Reynolds-stress anisotropy bij and 
of dissipation-rate anisotropy b^.^ ) , demonstrating how such r elatively compact models 
can be highly successful. On the other hand, the S ( Sugal[2004} ) model of the tcl- family, 
is built using a much greater number of reducible tensors, with an effort to have few 
non constant coefficients. Nonetheless, when projected onto the irreducible polynomial 
representation basis (Tab. [3]) all the coefficients are functions of the inv ariants (in the 
prese nt ca se of Reynolds-stress ani sotropy only). Finally the LRR (jLaunder et al 



19911 ) models include the 



1975f) . GL ('Gibso n fc Laundei|[l978l ). and SSG ([Speziale et al 
anisotropy of dissipation, ie model (j>ij — 2pe5e.^- (Lumley 1978 

There exist relatively fewer proposals for pressure diffusion d- ?-* fM.4D. The a priori 
TMC assessment of the weakly- inhomogeneous proposal (|4.19p of Lumlev ( 1978[ ). using 
the DNS-computed triple-velocity correlations shows that this closure cannot correctly 



model both d^J and d^fj , in wall turbulence, the later normal compon ent d'^fj bein; 



(p) 



I 



predicte d with the wrong sign. This issue notwithstanding, the models of IValletl (|200 
and Fu ([l99i) give the best predicti on of the DNS d ata f ij4.5p . Since both are based on 
some form of the proposal (j4.19p of Lumlev ( 19781 ). they require improvement by the 
addition of some other term in the representation of either p'u'^ or c?^^-* directly. 



The a priori assessment shows that all of the models for (^ H4.3p and for d^"^ (fj 
are unsatisfactory in the very-near-wall region (y"*" S 3) where they all return nearly 0- 
values, contrary to the DNS data. This problem does not necessarily affect the global RSM 
model performance, since in this region (y+ ^ 3) 4>ij and cancel one another, so that 



the prediction of 11^ — <pij + d^^^ (13.21) , which is the term appearing in Reynolds-stress 

which goes 



transport (13. ip , is not necessarily unsatisfactory, with the exception of 11 
to much slower than the other components. Away from the wall, where (jjij and d^^ 



do not cancel one another, there are reasons to believe that pressure diffusion plays 
an importa nt role in the predict i on of shear layer entrainment and should be modelled 
separately (jPemuren et aZ.lll996HValletl[2007l ). 

Future work will concentrate in extending the wall-echo DNS database to higher Reynolds 
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numbers, and in using the data to improve near-wall modelling of the pressure correla- 
tions. 
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Appendix. Green's functions 



As discussed in we have to solve the generic modified Helmholtz (iCheng et 
20061) equation 



(Al) 



with associated boundary-conditions at y G {-^-, L^}, for the complex-valued function q : 
R — !• C, Q{y) being a given complex-valued function Q : M — > C, and k £ M>o. Notice 
that there is no loss of generality in assuming k £ R>o because (lA ip depend s on and is 
there fore independent of the sign of k. The general method of solution (Ben der fc Orszae 
1978[) of the linear ode (|A ip is based on the determination of the appropriate Green's 



function G{y,Y; k), solution of 



9y2 



K'Giy,Y;K,)=d{y-Y) 



lim 

e->-0+ 



f^G r)G 

^{Y + e,Y;K)-^{Y-e,Y;^) 

G{y,Y;K)eG'>{L^,L+) 
G{y,Y;^)£G\L_,L+)\{Y} 



(A 2 a) 

(A 2 6) 

(A 2c) 
(A2rf) 



so that 



9(2/; K-) 



G{y,Y;K)Q{Y;K)dY 



(A3) 



satisfies (jA ip because of (jA 2p . The Green's function satisfies the symmetry (reciprocity) 
condition 



G(y,y;^) = G(y,y;At) Vy, F G (L_ , L+) 



(A 4) 



because the modified Helmholz oper ator \d'f {■) — k'^{-)] is a self-adjoint linear differential 
operator (|Inc&.1926:.Courant &: Hilbert.il953[ ). As shown in lBender &: Orszag, (1978.) . the 
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general solution of (|A 2[) is 

A^{Y-n)qi{y-K) + A2{Y-n)q2{y;K) y <Y 
q2{Y-K) 



G{y,Y-n) 



Ai{Y-k)- 
+ [A2{Y;n) + 



[W{qi,q2)]{Y-K 
[W{q,,q2)KY;^i) 




(A 5) 



y>Y 



where qi{y; k) and 92(2/; k) are 2 linearly independent solutions of the homogeneous equa- 
tion q"{y) - Kq{y) = 0, 



[W{qi,q2)]{y;K) :=det 



qiiy,!^) 92(2/; k) 



^0 



(A 6) 



is the Wronskian (•' in (jA 6p denotes differentiation by y, k being a parameter), which 
is 7^ iff the 2 solutions are linearly independent. The functions Ai{Y; k) and A2{Y; k) 
are determined by the boundary-conditions. The cases k ^ and k — are fundamen- 
tally different, not only because the 2 linearly independent solutions of the homogeneous 
equation q"{y) — K,q{y) — differ, but principally because k = is a characteristic 

number (eigenvalue) of (jA viz q"{y) — nq{y) = (l"{y) — with homoge neous Neu- 
mann bo undary-conditions admits the constant function as an eigensolution ( Incd 1926 : 



Courant fc Hilbert.lQSaV 



A.l. K^O 



When K 7^ we have 



gi(y;K^O) =e+«^ 
q2{y;n^{)) ^e-^y 

[W{qi,q2)\{y;n^Q) =det 
so that solution (jA 51) reads 



-Ke~'^y 



^2k ^ 



(A 7a) 
(A7&) 

(A 7c) 



Ai {Y; K)e+''y -f A2 ; K)e-'^y y <Y 



G{y, Y;k^O) = ^ Ai{Y; K)c+''y +^2(1"; K)e-''y 

+ — sinhfrefy — Y)] 

K 



(A 8) 



y>Y 



The functions Ai{Y;k) and A2(F;k) are determined so as to satisfy the boundary- 
conditions. 

A.1.1. GKin.{y,Y;K^O) 

Kim (Il989l) solves for the homogeneous Neumann boundary-conditions (j2.66p 



— (y — ± — .Y\k] 

dy 2 ' ' ' 







dy 2' ' 







Straightforward calculations yield the functions Ai{Y]n) and A2{Y;k) in 
satisfy (jA9ap . giving 



GKim(2/,r;K^0) = - 



cosh[K(ij^ - |y - Y\)] + cosh[K(y + Y)] 
2Ksinh K,L„ 



(A9a) 
which 

(A96) 
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This if-less expre ssion, which highhghts the symmetry condition (jA4p . is equivalent to 
the form in iKim (|1989,. (8), p. 440), written here using dimensional variables 



cosh[K{^ Ly — Y)]cosh[n{^ Ly + y)] 



K sinh kL,, 



cosh[K(iLy — y)] cosh[K(iLy + ¥)] 



K sinh kL^ 



y<Y 

y>Y 



(A 9 c) 



because using well-known identitie^ (jHarris fc Stockeil 1998 . p. 249) 

2 cosh[K(i_Ly — Y)] cosh[K{^Ly + y)] = cosh[K{Ly — Y + y)] + cosh[K(— F — y)] 

= cosh[K(L2^ -\y-Y\)]+ cosh[K{y + Y)] 
yy<Y {A9d) 
2 cosh[K(iLy — y)] cosh[K(^Lj, + Y)] = cosh[K(L.y — y + Y)] + cosh[K(— y — Y)] 

= cosh[K(Ly - |y - Y\)] + cosh[K(y + Y)] 
yy>Y (A9e) 

A. 1.2. Accommodating the boundary- conditions (k ^ 0) 

To take into account normal gradients dyq{y = ±^Ly;K ^ 0) = B+(k ^ 0) ? 7 a-t 
the channel walls, corresponding to the Stokes pressure (12 3[ 12 61). IChang et al\ ( 19991) 
superpose a boundary-condition q^c term to the solution of KimI (|l989 ) 



q(2/;«;^0,B±) - 

d\{y-K + (\B±) 
dy^ 



GKim(y, AC ^ 0)g(r; k ^ Q)dY + q^^{y- k ^ 0, B±) (A 10a) 
K'q{y;K^{),B±)=Q{y;K^Q) y e i-^Ly,+^Ly) (A 106) 



d^qBc{y;K^O,B±) 2 / /nn \ n ^^ir ,ir 



dy 



) (A 10c) 



^(y = ±iL,; ^ ^ 0, Si) ™ ^(y = ±iL,; At ^ 0, i?±) = S±(^ ^ 0) (A lOd) 

By (lA 3[ lA 9[) . the superposition (lA lOaj) is the solution to ()A 10 6p . with boundary- 
conditions (lA 10(ip . The solutio n of (lA 10 cl) with b oundary-conditions ()A 10d|) is readily 
obtained by direct integration ( Chang et al. Il999l) . and reads 



gBc(y;K 7^ 0,B±) = 



5+ cosh[K{^Ly + y)] — B_ cosh[K{^Ly — y)] 



K sinh kL 



(All) 



A. 1.3. Freespace Green's function G<xs{y, F; k ^ 0) 

The freespace Green's function is directly obtained from the general solution (jA 81) by 
requiring that the solution should tend to as y — > ±oo, which is equivalent to 



lim K 7^ 0) = 



lim Gffl(y,r;Ac^O) -0 VF G M \ {y} (A 12) 

>±oo 



Using (|A 12p to determine the functions Ai(Y] n) and A2(Y] k) in (lASp readily gives 



G2,(y,r;K^0) = 



2k 



t cosh(a -\-h) = cosh a cosh 6 -I- sinh a sinh 6 cosh(a - 



(A 13) 

cosh o cosh b — sinh o sinh b 
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A. 1.4. Upper/lower wall Green's Junctions G±{y,Y; k ^ 0) 

These Green's functions correspond to the problems where the domain of interest 
is (— oo, ^Ly) (isolated upper-wall influence; G+) or {—^Ly, +oo) (isolated lower-wall 
influence; G-), ie they are halfspace problems (±y e (— oo, \Ly)). For each of these 
problems only one of the walls is present. The associated boundary-conditions are 0- 
gradient at the wall of the problem, and bounded influence at the infinity of the problem, 
ie in terms of Green's functions 

G+(y,y;K 7^ 0) satisfies (|A2]) <^ dMI (A 14a) 

(A 146) 



dG 



lim G+iy,Y;K^O) =0 VF e M \ {y} (A 14c) 

to study the virtual influence of the upper wall, and 

G_(y,y;K ^ 0) satisfies (IX2| <^ (|X3l) (A 15a) 
dG 

-^{y^~\Ly,Y-K^Q) = Q (A156) 

lim G_(y,y;K 7^ 0) = VreM\{y} (A 15c) 

to study the virtual influence of the lower wall. By (jA 3| . boundary-conditions (|A 14 &|, 
lA 15fe|l imply 0-gradient at the wall of the problem (upper for G+ and lower for G_). 
Boundary-conditions (jA 14cl I A 15 cP ensure boundedness at infinity for the halfspace cor- 
responding to each problem. Straightforward computation of the functions AiiY; k) and 
A2{Y; k) in (|A8p yields 

G_(y,y;^^0) = y,Y e[-^Ly,+^) (A 16a) 

G+iy,Y;K^O)^ y, F e (-cx), +iL,] (A 166) 

which by (lA 13p reads 

„-K{Ly+y+Y) 

G_(2;,r;«^0) = G2j(y,r;«^0) y,Y e [-^Ly,+^) (A 17a) 

ZKi 

G+{y,Y;K^O)^G^{y,Y;K^O)-- ^ y^Ye{-^,+^Ly] (A 176) 

ie the presence of a single wall in any of the 2 halfspace problems induces an additive 
correction to the freespace Green's function, the wall-echo Gw± {y, Y; k ^ 0) 

G^^{y,Y;K^O) :^ G±{y,Y; k ^ 0) - G<s{y,Y; k ^ 0) (A 18) 



A. 1.5. Accommodating the boundary- conditions for the halfspace problems (k ^ 0) 

To take into account normal gradients dyq{y — ±^Ly;K) = B±{k) 7^ at the wall 
for each of the halfspace problems of i 3A.1.4l we may proceed exactly as in i jA.1.2| by 
superposing a boundary-condition qBc± term to the corresponding solution (|A 15[) with 
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G^{y, Y;k^ 0)Q{Y; k ^ {))dY + q^^,_ {y; k ^ 0, S_) (A 19a) 



^^^^^^^^^^^i^-«25(2/;^^0,B-)^g(y;^^0) ye(-iL„+oo) (A 196) 



K^qBc.iy;K,B-) =0 y e {-^Ly,+(X)) (A 19c) 



^Ly;K^O,B_)=B_{K^O) (Al9d) 



lim gBc„(2/;'« ^ 0,B_) = (A 19e) 



to study the virtual influence of the lower wall boundary-condition on the solution, and 



q{y;K^O,B+) - 

d^qiy;K^O,B+) 
dy"^ 



G+{y, Y-K^ 0)Q(r; k ^ Q)dY + Qbc^ {y- ^ 0, B+) (A 20a) 



-n\{y-K^O,B+)=Q{y-K^O) ye{~^,+\Ly) (A 206) 



d^qBc+iy;K^ O-.B^ 
dy"^ 



n'qBcM,K^Q,B+)^0 y^{-^,+\Ly) (A 20c) 



g(y = +iL,; « ^ 0, i?+) ™ ^ +iL,; ^ 0, i?+) = ^ 0) (A 20d) 

lim (7bc+(2/;a«^0,B+) = (A20e) 

to study the virtual influence of the upper wall boundary-condition on the solution. By 
(|A Isp . the superposition (|A 19 aP is the solution to (jA 196p . with boundary-conditions 
(jA 194 IAl9ep . By (|A14p . the superposition ()A20aP is the solution to (|A206p . with 
boundary-conditions (|A 20 d I A 20 ep . Straightforward computation yields 



gBc± (y; 0, B±) = ±rie-''(2^«T^) 



A.2. K = 

When K = (jA ip becomes a Poisson equation (jKatz fc PlotkinlHool 

d^q{y;K = 0) 



(A21) 



dy2 



(A 22) 



The limits as k — > of the solutions obtained in i jA.ll for k ^ are all singular, as can be 
easily verified by straightforward coi nputation contr ary to the case of the same problem 
with Dirichlet boundary-conditions ( Zauderen 2006[ ) . This is related to the existence of 
solutions of (jA22p with Neumann boundary-conditions at y — ±^Ly, because integrating 
((X22|l yields 



|(+lL,;. = 0)-|(-iL,;. = 0) = 



Q{Y;k=^ 0)dY 



(A 23) 



implying that (jA 22p with gradient boundary-conditions q'{±^Ly, k = 0) — B± can 
only be solved iff the compatibility condition (jA 23P is satisfied by the integral of the 
source-term and the boundary-conditions. 
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In the particular case of homogeneous Neumann boundary-conditions B± = 0, which 
is the most important since B± ^ is only concerned with Stokes pressure p. 31 12.6p . 
the homogeneous equivalent of equation (jA 22[) . q"{y,K = 0) = 0, admit s any con- 
stant function as soluti o n, ie any constant functio n is an eigenfunction (jlncd Il926 : 
Courant fc Hilbertlll953t iMvint-U fc DebnathI [2007h of the homogeneous equivalent of 



(|A 22p . q"{y, k = 0) = 0, with boundary-conditions q'{±^Ly, k = 0) = 0, and k = is a 
characteristic number of the homogeneous equivalent of (jA ip . q"{y,K) — K^q{y,K) = 
with boundary-conditions q'{zt^Ly, k) — 0. Therefore a solution only exist s provided the 
corresponding compatibility relation (Ince 1926 : Courant fc Hilbertlll953l) holds, ie the 
integral in (jA 23p is equal to 0. This is analogous to the 3-D compatibilit y relation for the 
Poiss on equation with Neumann boundary-conditions ( Ockendon et a/., ,2003; .Zaudered 



A.2.1. GKnn{y,Y;K = 0) 

In this case (j2.3M \2.Qb[ 0-gradient boundary-conditions apply on both walls (y 



Assuming that the compatibility relation 



(A 24) 



holds, there are 2 different ways for constructing a modified Green's function (* Coura nt fc Hilbert 
[l953.;,Zaudcrcr 200l| for solvi ng (lA 22|) w ith homogeneous Neumann boundary-conditions. 
The first method described in lCourant fc HilbertI (jl953l ) makes explicit use of the eigenso- 
lution of (jA 22j) in a modified definition of the problem (jA 2j) for calculating the modified 
Green's function whi ch satisfies the boundary-conditions, but has also to modify (jASjK 
The s econd method (jZaudereii 120061 ) which leads to the Green's function used bv iKim 
( 1989[ ) uses the standard definition of the Green's function (jA 2\ along with the stan- 
dard integral solution (jA 3p . In this latter case however the boundary-conditions are not 
included in the solution but are obtained because of the compatibility relation (jA 24^ . 

Hence assuming that the compatibility relation (lA 24^ holds, so that (jA 9aP for k = 
may be satisfied, the Green's function corresponding to (jA 22p satisfies 



dy"^ 



5{y^Y) 



lim 

£-i-0+ 



Kim 



dy 



dGv 



dy 



-(r-e,y;K = 0) 



GKim(y, Y-n = {))(^ C^'i-^Ly, +\Ly) 
GkUv^Y- k) e C\-\Ly, +\Ly) \ {Y) 



(A 25a) 

(A 256) 

(A 25c) 
(A25d) 



It is therefore made up by 2 straight lines joined together at Y , since by (jA 25 ap dyyG{y, Y; k 
0) = 0\fy ^Y. If we further admit that GKim{y, F; k = 0) = GKim{Y, y; k ^ 0) (iXIl) 
as expected because of the self-adjointed ness of the the opera tor [rfyy(-)] dncd Il926t 
ICourant fc HilbertI [l953l ) . we must have (jKatz fc PlotkinI Il99l[ ) GKim{y,Y;K = 0) = 
GKimdy - F|; K = 0) so that 



GKim(y,'5^; k = 0) = 



Y\ 



(A 26) 



The factor ^ is required to satisfy (|A25&p . and c € M is a constant. The value of c has no 
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influence whatsoever on the solution (|A 3|) because of the compatibility condition (|A 24|) 
The simplest choice c = corresponding to the reasonable situation GKim(y, y, n — 0) = ^ 
is made 

GKim{y,Y;^^0)^^y^ (A 27) 

and we have 

^(±iL,; . = 0) f^^' ^{y, Y;. = 0) Q{Y- . = 0)dY 

r^"'QiY;.^0)dY^O (A 28) 



Notice that in this case it is the compatibility relation (jA 24[) . and not the Green's function 
alone which is responsible for satisfying the boundary-conditions. 

A. 2. 2. Accommodating the boundary-conditions (k — 0) 

Finally, the solution of the boundary-conditions problem (required for P(t-)) can be 
obtained directly by integrating 

ABc(.;. = 0,i^.,i?-) ^^ ..HX„+W (A 29.) 

= ±ii„;K = 0,B+,B) = B±(k = 0) (A 296) 

dy ^ 

The solution exists iff 

= 0) = = 0) (A30) 

and is obviously a straight line, defined up to an additive constant 

<ZBc(y;K = 0,B_,B+) = iS-y+iS+y + ^BCo (A31) 

where by (|A 3ip (7bc(2/ = 0; k = 0) = q^ca is the solution at centerline. 

A. 2. 3. Freespace Green's function G^{y,Y; k = 0) 

The solution for the freespace Green's function G<i3{y,Y; k = 0) is obtained exactly in 
the same way as for GKim{y, Y; k — 0) f! jA.2.1[) . yielding by (jA 27[) 

G<x,{y, Y;k = 0) = GKnniy, Y;k = 0) = (A 32) 

Again, the sources Q{y\n = 0) must satisfy the compatibility condition (jA 24[) to ensure 
boundedness as |?;| — )• oo 

lim / G<a{y,YK^Q) Q{Y\K^Q)dY 
^ lim l^y"^' "g(y;K = 0)dy^ M| VreM\{y} (A33) 

t / ( ^ +cj Q{Y-n^ Q)dY = / ^ \y—AQ(Y-^, = Q)dY+cl k = 0)dF 



IIA32t 
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which is iff (|A24l) holds. 



A. 2. 4. Upper /lower wall Green's functions G±{y,Y; k — 0) 

The solution for G±{y,Y; n — 0) is again obtained in exactly the same way as for 
GKnniy,Y; ^ = 0) (ElID, yielding by ^EM 

G± iy, y ; « - 0) - G^T (y, K = 0) = GkUv, Y;k^O)^ (A 34) 

The boundary-conditions for the gradient at the wall y — ^^Ly, and for boundedness 
at ±cx), are again satisfied by relations analogous to (|A28[ lA 33p . iff the compatibility 
relation (|A24p holds. 



A. 2. 5. Accommodating the boundary- conditions for the halfspace problems (k ~ 0) 

Finally, the solution of the boundary-conditions (required for P(^)), in the case of the 
halfspace problems, can be obtained directly by integrating 

Ty e (-iLj,,+oo) (A 35a) 

0,B±) = B±{k = 0) (A 356) 

The solution is obviously a straight line of slope B±{k — 0), and is bounded as \y\ — ?> oo 
iff the compatibility condition 

B±{K^O)=0=^qnc^{y;K = 0,B±)=0 (A 36) 

holds. 



dy2 

dqBc± 
dv 



(y 



A. 3. Approximation error in the method of images 

One interesting application of the halfspace proble ms is that, by superp osition, they 
give the solution obtained by the method of images (jManceau et aLl 120011) . Indeed the 
method of images corresponds to applying the freespace Green's function on the domain 
y G {—^Ly,+^Ly), and on 2 ghost domains y G {—^Ly,—^Ly) U {+^Ly,+^Ly), with 
appropriately reflected source-terms 

QMWL{y;n) =Q{-Ly - y;K) y e i-^Ly,-^Ly) (A 37a) 

QMWL(y;K)= <9(y;K) y ^ i-^Ly^+^Ly) (A376) 

QuwLiyiK) ^Q{+Ly - y; k) y £ {+lLy,+^Ly) (A 37c) 

and th e approximation to the solution obtained by the method of images used bv lManceau et cd. 
(|200lHt l is 

quwhiy^ii)^ / G9j(y,y; k) Qmwl(>'; k) (A38) 

It is straightforward to construct a kernel GMWL{y,Y;K) which applied on the source- 
terms in the actual computational domain y G {—■^Ly, —^Ly) will return the method of 

f Altliougli'Manceau et al] (|200lf ) neglect qBciy, f^, B±), this can be added to quwhiv, n) (|A 38|l 
to obtain tlie approximation of the complete solution (|A 10a[) . 
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images, because, defining Y± ±ij, — Y, we have 



Quwhiy; k) - QBciy; k,B±) = / G<s{y, Y; k) QuwhiY; k) dY 

2 



l lA37l .l lA38l 



-i r 
2 



G^^{y,Y;>i) Q{~Ly-Y-K) dY 
G^{y,Y;K) Q{Y;k) dY 
G^{y,Y-K)Q{+Ly-Y;K) dY 
G<a{y, -Ly -Y^;k) Q{Y^; k) dY^ 
G<aiy,Y;K) Q{Y;k) dY 
G^{y,Ly-Y+-K)Q{Y+;K)dY+ (A 39) 



yielding 



quwdy, = I Guw^iy,Y;K,) Q{Y;k) dY (A 40a) 

Gmwl(2/, Y; k) ■.^G<n{y, Y; k) + G^xi{y, -Ly ~Y;k) + G<n{y, Ly - Y; k) (A 406) 

and by ((XT3)) 

fXr^l p-i^\y-Y\ (.-ii\Ly+v+Y\ -K\Ly-y-Y\ 

GuUy,Y;K^ 0) ^ - - £ _ 2 _ (A 40c) 

Notice that (jA 40p give^ the solution by integrating the sources in the actual domain 
Y S [—^Ly,—^Ly] with an approximate Green's functions GMw-L{y,Y;K) (|A406[) . This 
representation (|A 40[) readily yields the approximation error of the method of images, by 
comparison with the exact solution (jA 96p obtained using GKini{y,Y', k). If we consider 
the actual domain y,Y E [—^Ly, +^Ly], we have 

-^Ly<y<+lLy \ ( Ly + {y + Y)>0 
-\Ly<Y<+\Ly ]=^\Ly-{y + Y)>Q ^^'^^"^ 



t l|X40)) is the 1-D equivalent of (|Manceau et aZ.ll2001l . (4.2), p. 314) 
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2k 

Q-K{Ly+y+Y-^ 

2k 



Vy,y e [-iL^,+iLy] (A416) 



(A 41c) 



+ G-{y,Y;K^O)-G^iy,Y;K^O) 

G^_iy,Y;K^O^KM) 

+ G+iy,Y;K^O)-G<siy,Y;K^O) 

^ ^ ~ ' 

G„+(y,y;K^O) ^KM 

yy,Y e[-^Ly,+^Ly] 

ie, the images problem (jManceau et al. I l200lh is equivalent to adding to the freespace 
Green's function an independent wall-echo correction for each wall, the wall-echo correc- 
tion for the upper (lower) wall being exactly the wall-echo of the halfspace problem with 
only the upper (lower) wall present fi iA.1.4[) . Notice that by (|A34[) Gu,_ iy,Y; k = 0) = 
Gw+ {y, F; K = 0) = so that (IA41c|) holds Vk G M>o. 

By (jA 961 IA41bp we may directly compute the error made when approximating the 
Green's function by the method of images 

^ (G„w.(y, Y;n^O)- G^Uv. - ^ 0)) 

cosh K\y-Y\^ (f~>^Lv cosh K(y + Y) 



2KLy sinh nLy 



(A 42) 



\fy,Y e[-^Ly,+^Ly] 



where (jA 421) was made nondimensional by dividing the Green's functions by Ly. Obvi- 
ously 



y,Y e[-^Ly,+^Ly] 



0<\y~Y\<Ly 
0<\y + Y\<Ly 



< coshKly — Y\< cosh nLy 
< cosh + y I < cosh nLy 

(A 43) 

since coshx = cosh|a:| is positive ans strictly increasing V|a;| G M>o. HenccQ 



— max 



GMWLiy,Y;K)~GKim{y,Y;K)j < 
2KLy sinh nLy 



(A 44) 



EuWhil^Ly 



showing that the approximation error made by the method of images is a function of the 
nondimensional (outer scaling) wavenumber KLy (Fig. 2]). 



t notice also that by (|A42ll GymL{y,Y; k ^ 0) -GKim(2/, F; k / 0) > Vy.F £ [-^Ly,+^Ly] 
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